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ABSTRACT 


ON  THE  PHYSICAL  REALIZABILITY  OF  BROAD-BAND  EQUIVALENT 
CIRCUITS  FOR  WIRE  LOOP  AND  DIPOLE  ANTENNAS 


Results  of  a  study  of  the  feasibility  of  using  Singularity  Expansion  Method 
(SEM)  data  to  synthesize  broad-band  equivalent  circuits  for  loop  suid  wire 
structures  are  reported.  The  positive  realness  of  admittances  derived  are 
studied,  on  the  basis  of  both  pole-pair  groupings  and  eigenmode  groups  within 
the  SEM  formalism.  Synthesized  circuits  are  given  for  a  center-driven  and  a 
quarter-driven  straight  wire  antenna,  and  for  a  circular  loop.  Results  are 
compared  with  those  derived  by  means  of  numerical  solution  of  the  time- 
dependent  integral  equations  for  each  structure.  The  sensitivities  of  the 
response  to  changes  in  the  values  of  poles  and  to  circuit  element  valioes  are 
considered. 
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CHAPTER  I 


INTRODUCTION 

1.1  Description  of  SEM 

(fuch  attention  has  been  focused  on  the  use  of  the  singularity  expan¬ 
sion  method  (SEM)  to  solve  broad-band  transient  electromagnetic  field  problems 
Involving  scatterers  since  Its  formalism  was  developed  by  Baum  [1,2]  In  1971. 
SEM  formalism  grew  from  the  experimental  observation  that  the  transient 
response  of  scatterers  Is  dominantly  composed  of  one  or  more  damped  sinusoids. 
Since  such  damped  sinusoids  correspond  to  conjugate  pole-pairs  In  the  complex 
frequency  plane  (or  s-plane) ,  Baum  postulated  that  a  complete  representation 
of  the  transient  response  could  be  obtained  by  knowledge  of  all  s-plane  sin¬ 
gularities  of  the  object's  response.  Such  s-plane  singularities  potentially 
Include  poles,  essential  singularities,  branch  cuts,  and  entire  functions.  A 
complete  SEM  representation  of  a  scatterer  consists  of  a  set  of  singularities, 
the  modal  current  distribution  associated  with  a  singularity,  and  a  normaliz¬ 
ing  coefficient.  Marin  [3]  showed  that  perfectly  conducting  objects  In  free 
space,  such  as  are  dealt  with  In  the  present  work,  have  only  pole  singular¬ 
ities.  Hence  only  pole  singularities  are  considered  In  this  work. 

Because  the  SEM  approach  reduces  the  problem  of  finding  a  transient 
response  for  a  given  object  to  that  of  finding  the  object's  singularities,  the 
associated  modal  distributions,  and  the  coupling  factors  which  welgiht  a  given 
singularity's  contribution  to  the  response  for  a  given  exciting  waveform,  we 
see  chat  SEM  offers  a  more  compact  and  efficient  way  of  representing  transient 
responses  than  previous  methods.  These  methods  Include  numerical  solutions  of 
dlfferentlal/lntegral  equations  In  the  time  or  frequency  domain,  and  Integral 
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operator  taodal  approaches.  The  connection  of  the  latter  approach,  sometimes 
referred  to  as  the  eigenmode  expansion  method  (EEM) ,  to  SEM  has  been  demon¬ 
strated  by  Baum  [4].  In  the  present  work  extensive  use  of  EEH  in  conjunction 
with  SEM  is  made.  The  utility  of  SEM  lies  in  the  low  and  intermediate  fre¬ 
quency  ranges,  and  is  particularly  useful  for  studying  the  response  of  scat- 
terers  to  the  nuclear  electromagnetic  pulse  (EMP)  because  of  the  presence  of 
frequencies  whose  wavelengths  range  from  long  to  short  compared  to  scatterer 
dimensions  in  such  pulses.  In  addition  to  compactness  and  efficiency,  SEM 
offers  greater  physical  insight  into  a  problem  than  other  methods. 

1.2  Scope  of  Present  Work 

The  motivation  for  the  present  work  derives  from  the  insight  that 
it  may  be  possible  to  construct  passive  RLC  (resistor.  Inductor,  capacitor) 
circuits  from  the  pole  and  residue  data  of  SEM.  Such  circuits  might  then  be 
used  to  predict  the  transient  current  response  of  an  object  to  various  wave¬ 
forms.  This  idea  has  been  formalized  by  Baum  [5,6]  and  the  present  work 
centers  on  the  feasibility  of  this  approach  for  the  construction  of  some 
realizable  equivalent  circuits  for  radiation-gathering  structures  -  in  par¬ 
ticular  on  the  positive  real  considerations  of  realizability  for  convenient 
formalisms. 

In  a  previous  work  Hess  [17]  has  performed  an  EMP  coupling  analysis 
using  transfer  Impedances  derived  from  SEM.  Schaubert  [16]  has  recently  con¬ 
structed  lumped-element  equivalent  circuits  for  a  center-driven  dipole  and  a 
Yagi  antenna,  using  experimentally  derived  SEM  data.  Schaubert  extracted  SEM 
data  by  means  of  Prony's  algorithm  from  the  antenna  terminal  voltage  wave¬ 
form  due  to  a  step-like  excitation,  and  used  this  data  to  construct  the  total 
impedance  of  the  antennas.  Schaubert  then  used  Brune's  synthesis  method  to 
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derive  equivalent  circuits.  The  present  work  uses  theoretically  derived  SEM 
data  and  addresses  the  problem  of  realizing  the  antenna  Impedance  by  Indi¬ 
vidual,  transformerless,  circuit  modules  which  can  be  related  to  the  SEM  data. 

Baum  [5,6]  suggested  the  construction  of  equivalent  circuits  on  a 
pole  at  a  time  basis,  a  conjugate  pole-pair  basis,  and  an  eigenmode  basis,  and 
gave  generalized  formal  circuits  for  each  type  of  synthesis.  The  realizabil¬ 
ity  of  the  formal  circuits  was  not  treated  In  detail.  The  formal  development 
of  SEM  leading  to  equivalent  circuits  Is  left  to  Chapter  II.  The  case  of 
conjugate  pole-pair  synthesis  Is  examined  briefly  and  related  to  earlier  work 
by  GulUemln  [7],  In  order  to  Illustrate  some  of  the  salient  aspects  of  the 
current  work.  A  more  detailed  consideration  of  this  and  related  problems  Is 
given  In  Chapter  III. 

The  physical  realizability  of  a  given  Impedance  or  admittance 
quality  hinges  on  whether  or  not  that  quantity  Is  a  positive  real  (PR)  func¬ 
tion  of  the  complex  frequency  s.  A  positive  real  function  Is  defined  as  a 
function  that  Is  real  for  s  real.  Is  analytic  In  the  right  half-plane,  and 
whose  real  part  Is  positive  along  the  Jcu  axis.  Additionally,  any  poles  on 
the  Ju  axis  are  simple  with  real  and  positive  residues.  It  Is  Instructive  to 
study  the  conditions  under  which  a  conjugate  pole-pair  exhibits  positive  real¬ 
ness  ("FRness’’) . 

The  admittance  associated  with  the  n^^  SEM  pole  Is  written  as 

a 

■  .  (1.1) 

n  s  - 

n 

where  ■  o  +  Ju  Is  the  complex  pole  and  a  “a  +  JS  Is  the  complex 
nn-'n  rr  nn'^n 

residue.  A  modified  form  of  this  admittance,  which  has  the  property  of  being 
zero  at  zero  frequency.  Is  written  as 


so  that  Y^(o)  “  0  is  physically  consistent  with  the  terminal  admittance  of  a 
simply-connected  scatterer. 

SEM  poles  are  either  purely  real  or  occur  In  conjugate  pairs.  The 
conjugate  pair  case  Is  the  more  general  one.  The  admittance  of  the  n^^  con¬ 
jugate  pole-pair  Is  then.  In  tinmodifled  form. 


Y=P(s) 


(1.3) 


or 


2  Real(a^)s  -  2  RealCa^s 

(s)  -  —5 - 2 - ^  .  (1.4) 

s  -  2  Real(s  )s  +  si 
n  n' 


In  this  form  we  see  that  the  coefficients  of  Y^^(s)  are  real  valued,  as  they 
must  be  for  circuit  realizability.  The  other  criterion  necessary  to  ensure 
realizability  of  Y^ (s)  is  the  non-negatlveness  of  its  real  part  along  the 
ju  axis. 

Gulllemln  [7]  showed  that  for  an  unmodified  admittance  function 


formed  by  conjugate  pole-pairs  to  be  PR,  certain  conditions  must  be  met.  If 
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d,  ■  2o_ 


b,  -  -2o 


-2  Real(aV) 
n  n 


I  |2  2.2 

S  ■  O  +  (1) 

'  n'  n  n 


-2(a  o_  +  6_u  ) 
n  n  n  n 


(1.6) 


We  see  that  and  b^  are  positive  for  a  pole  in  the  left  half-plane,  while 
and  d^  may  be  of  either  sign.  In  order  to  find  what  conditions  need  to 
apply  to  dj^  and  dQ  for  Y^^(s)  to  be  realizable,  we  need  to  find  its  real  part, 
and  check  for  non-negatlveness  at  s  ••  Ju.  This  leads  to  the  condition 

dgbQ  +  (djbj^  -  >  0  .  (1.7) 

This  requires 

dj^bj^  >  dQ  >  0  ,  (1.8) 

which  leads  to  the  PR  conditions 

1 

-a  a  +  g  u  >0  (1.9) 

ti  ti  II  n 

I 

and 


-a  o  -6m  >  0  . 

n  n  n  n 


(1.10) 


We  can  combine  these  conditions  Into  a  single  expression  which  reflects  the 
necessary  restriction  on  the  residue  for  a  PR  function  to  result.  This 
expression  Is 

(1.11) 


This  condition  Is  Illustrated  graphically  In  Figure  1.  These  and  similar  PR 
conditions  are  used  extensively  in  Chapter  III.  For  Y^^(s)  to  be  PR,  the 
residue  must  lie  within  the  shaded  portion  of  the  a-plane.  Unfortunately,  It 
Is  seen  that  most  SEM  poles  and  residues  fall  to  meet  this  criterion,  and 
ways  of  circumventing  this  difficulty  are  devised.  The  parameter  tests  which 
are  required  for  modified  conjugate  pole-pairs  have  been  derived,  and  are 
presented  In  Chapter  Ill. 

In  Chapter  II,  the  formal  development  of  SEM  in  conjunction  with 
EEM  as  applied  to  equivalent  circuits  Is  given,  and  the  sources  of  SEM  data 
used  In  this  study  explained.  Chapter  III  deals  with  PR  considerations  of 
ftoictlons  derived  from  SEM  data  on  conjugate  pole-pair  and  eigenmode  bases, 
and  examines  pole/resldue  error  considerations.  Chapter  IV  encompasses  some 
synthesized  circuits  for  wire  and  loop  objects,  and  compares  the  results 
obtained  from  these  circuits  and  analyzed  by  a  SCEPTRE  circuit  analysis  pro¬ 
gram  to  results  obtained  from  a  thin  wire-time  domain  (TWTD)  program.  Also 
In  Chapter  IV  the  sensitivity  of  these  circuits  Is  considered.  Chapter  V 
gives  the  conclusions  reached  by  this  study,  and  points  out  areas  which  need 


further  consideration 


/ 


CHAPTER  II 

FORMAL  DEVELOPMENT  OF  EQUIVALENT  CIRCUITS 

2.1  Introduction 

In  this  section  we  consider  the  problem  of  the  broadband  admittance 
of  an  arbitrary  scatterer  or  antenna  at  a  gap  region  In  the  object.  The 
approach  Is  to  cast  the  object's  terminal  pair  formed  by  the  gap  Into  the 
form  of  a  Norton  equivalent  circuit.  Then,  using  SEM,  the  equivalent  admit¬ 
tance  Is  found  In  the  form  of  a  residue  series  In  the  complex  frequency 
variable.  Subsequently,  network  synthesis  techniques  are  applied  and  equiva¬ 
lent  circuits  constructed.  Three  different  generic  circuit  configurations 
are  considered,  corresponding  to  three  different  levels  of  grouping  the  SEM 
poles.  In  the  first  circuit,  the  admittance  formed  from  the  Individual  poles 
Is  the  basic  module.  In  the  second  the  pole-pair  admittance  formed  by  group¬ 
ing  a  conjugate  pole-pair  together  Is  the  basic  module.  Finally,  the  admit¬ 
tance  formed  by  grouping  the  poles  according  to  eigenmodes  with  which  they 
are  associated,  as  directed  by  EEM  results,  provides  the  basic  module.  The 
admittance  formed  by  eigenmode  grouping  we  term  a  "terminal  elgenadmlttance". 

In  the  last  part  of  the  chapter  the  sources  of  the  SEM  data  used  in 
this  study  are  discussed. 

2.2  Formulation  of  the  Norton  Equivalent  Problem 

Figure  2  gives  the  steps  required  In  defining  the  Norton  eqtilvalent 
for  an  active  circuit  with  a  single  port.  The  active  circuit  can  be  replaced 
by  a  single  equivalent  admittance  and  a  current  generator  In  parallel,  as  In 
Figure  2a.  In  Figure  2b  the  methods  for  finding  the  current  generator  and 
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ACTIVE 

CIRCUIT 
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Isc(  s) 


Figure  2q 
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CIRCUIT 


Isct*  ^ 


ACTIVE 
CIRCUIT 
SOURCES =0 


Figure  2b 


Figure  2:  Defining  the  Norton  equivalence  admittance  and 
short  circuit  current. 
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equivalent  admittance,  both  necessary  for  the  replication  of  the  active  cir¬ 
cuit  as  seen  by  the  port,  are  Illustrated.  The  short-circuit  current  is  found 
by  shorting  the  port  of  the  active  circuit,  and  the  equivalent  admittance  is 
found  by  setting  the  sources  of  the  active  circuit  to  zero  and  measuring  the 
admittance  seen  by  the.  port.  The  following  development  is  essentially  that  of 
Baum  [5].  It  is  reproduced  here  for  the  s.ake  of  completeness. 

The  single  port  scatterer  problem  is  defined  in  Figure  3.  In  this 

figure,  S  is  the  surface  of  the  scatterer,  all  of  which  is  not  shown.  The 
& 

gap  which  is  formed  by  parting  the  scatterer  with  a  plane  is  denoted  bv  S  . 

g 

a  is  a  unit  vector  defining  gap  orientation  and  A  is  the  gap  width, 

g 

Using  this  geometry,  we  desire  to  find  the  Norton  equivalent  repre¬ 
sentation  of  the  antenna  as  seen  by  the  gap  (feed-point).  The  method  is 
illustrated  in  Figure  4.  The  equivalent  admittance  is  found  with  the  aid  of 
Figure  4a.  The  first  step  is  to  specify  an  electric  field  at  the  feed  point. 
The  ratio  of  gap  current  response  to  the  voltage  associated  with  the  specified 
field  yields  the  terminal  admittance.  We  choose  a  field  that  is  divergence¬ 
less  in  the  gap  volume,  i.e.,  a  quasi-static  electric  field.  The  field  is 
represented  then  by 


E  (r,t)  =  +V(t)e  (r)  , 


(2.1) 


where  •  V(t)  is  the  gap  voltage. 

Having  defined  the  electric  field  across  the  gap,  we  need  a  suitable 
definition  of  the  current  I(t)  through  the  gap.  For  t)iis  we  integrate  the 
current  density  over  the  surface  S  of  the  gap  and  average  over  the  longitu- 
dlnal  direction  of  the  gap. 


Figure  4a 


Figure  4:  The  scatterer/antenna  admittance  and 
circuit  current  formulations. 
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I(t)  “  j  J(r,t)  •  e  (r)  ds  dz 

''s  » 

g 

-  <J(r.t)  ;  e  (r)>  (2.2) 

8  g 

« 

where  J(r,t)  is  the  surface  current  density,  z'*’’  denote  the  edges  of  the 

gap,  and  the  notation  <>  denotes  spatial  integration  over  the  cylindrical 

g 

surface  forming  the  gap. 

The  equivalent  admittance  of  the  gap  is  then 

(a)  -  <J(^,s)  ;  e^(^)>  (2.3) 

*  V(s)  V(s)  ®  8 

where  -  represents  the  two-sided  Laplace  transform. 

The  procedure  for  finding  the  short-circuit  current  is  detailed  in 
Figure  4b.  Here  the  sources  are  represented  by  an  elertric  field  incident  on 
the  antenna.  The  short-circuit  current  through  the  shorted  gap  is  then  cal¬ 
culated  from  Equation  (2.2),  with  E.  (r,s)  *  e  (r)V^  (s). 

me  g  inc 

The  present  work  is  concerned  entirely  with  the  equivalent  admit¬ 
tance  depicted  in  Figure  Aa  and  the  positive  real  considerations  of  that 
solution.  Therefore  the  derivation  of  the  short-circuit  current  generator  is 
not  given  here.  Readers  interested  in  this  subject  are  referred  to  Baum  [5]. 

2.3  Equivalent  Admittance  Circuits 

We  construct  formal  equivalent  admittance  circuits  using  the  singu¬ 
larity  expansion  of  the  current  j(r,s)  as  a  tool.  Three  such  circuits  are 
constructed.  These  circuits,  in  order  of  the  complexity  of  the  modules  from 
which  they  are  constructed,  are  composed  of 

) 


I 


/ 
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-  pole  at  a  tine  admittances 

-  pole-pair  admittances 

-  terminal  elgenadmlttances 

The  singularity  expansion  of  the  current  In  the  frequency  domain  Is 


[5] 


?(r,s) 


V(s)  j^Z 


Wnl^^^  - 


.-1 


•nl' 


J®( 


r,s)j  , 


(2.4) 


where  V(s)  Is  the  voltage  across  the  gap,  Is  the  class  I  coupling  coeffi¬ 
cient,  and  natural  current  mode.  The  Index  n  Indexes  over  the 

eigenmode  sets  of  poles,  and  the  Index  1  Indexes  the  Individual  poles  In  an 

/ 

eigenmode  set.  J  (r,s)  Is  a  possible  entire  function  contribution.  Several 
researchers  [8]  have  conjectured  this  contribution  to  be  zero  when  the  sun^ 
matlons  are  ordered  over  elgenmodes  as  Indicated  In  (2.4).  The  entire 
function  contribution  Is  not  considered  In  this  work. 

Inserting  this  form  of  J(r,s)  In  Equation  (2,3),  we  obtain 

Y^(s)  -  <[j  I  (2.5) 


or 


Y  (s) 

g 


n  1 


-  Z  Z  Y^i(s)  +  Y®(s)  , 
n  1 


(2.6) 


(2.7) 


where 


“nl 


s„(r)> 

8  g 


(2.8) 


and  Y  (s)  Is  the  entire  function  contribution,  if  present.  The  term  a  .  is 

Ux 

termed  the  residue  associated  with  a  pole  s^^. 

Equation  (2.7)  is  now  used  to  construct  formal  circuit  diagrams.  The 

first  of  these  is  given  in  Figure  5,  in  which  the  modules  are  formed  from  one 

pole  at  a  time.  In  this  figure  the  port  representing  the  antenna  gap  is  on 

«<• 

the  right,  with  voltage  V(s)  across  it  and  current  I(s)  into  it.  The  modules 
are  placed  in  parallel  as  dictated  by  Equation  (2.7),  and  the  entire  function 
modules  are  also  shown. 

Many  objects,  such  as  dipole  antennas,  esdilbit  a  zero  admittance  at 
zero  frequency.  It  is  desirable  to  modify  the  above  admittance  forms  to 
exhibit  this  property.  We  do  this  by  subtracting  off  the  value  of  the  admit¬ 
tance  at  zero.  Writing  the  admittance  in  this  form  gives 


Y  (s)  -  I  I  a 
°  n  i 


hi 


[(,  -  +  rjl  .  V'  (.) 


Z  Z  y'  (s)  +  Y*  (s)  , 
n  1 


(2.9) 

(2.10) 


where  is  termed  the  modified  admittance  module  by  Baum.  This  form 

differs  from  the  unmodified  form  by  only  the  addition  of  a  constant  term. 

The  equivalent  circuit  construction  with  modified  admittances  would  take  the 
same  form  as  Figure  5.  We  note  that  the  modified  form  is  inappropriate  to 
the  pole  at  zero  for  multiply  connected  objects  because  the  subtractive 
factor  is  undefined. 

Since  in  general  both  the  residues  and  poles  in  Equation  (2.6)  are 


complex,  pole-at-a-tlme  modules  may  contain  unrealizable  elements.  This  is 
demonstrated  in  the  next  chapter.  As  a  step  toward  the  development  of  real¬ 
izable  circuits,  modules  formed  by  grouping  conjugate  pole-pairs  together  are 
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studied.  This  can  be  done  for  both  modified  and  unmodified  admittances.  In 
the  unmodified  case  we  have,  for  the  pole-pair  module. 


(8) 


(2.11) 


where 


<1<^> 


^nl 

S  —  8 


nl 


(2.12) 


and  ^jj^(8)  ^8  the  conjugate.  Figure  6  Illustrates  this  circuit  arrangement. 
In  this  figure  admittances  whose  poles  have  no  imaginary  part  are  designated 
by  \q(s). 

The  final  circuit  considered  consists  of  terminal  eigenadmlttance 
modules,  and  Is  represented  In  Figure  7.  Here  all  poles  belonging  to  an 
eigenmode  are  gathered  into  one  module,  and  the  complete  circuit  is  the 
parallel  sum  of  these  modules.  Again,  either  modified  or  unmodified  forms 
may  be  used.  The  module  for  the  n^^  unmodified  eigenadmlttance  is 


(2.13) 


\ 


In  Chapter  III  the  FBness  of  these  different  circuit  constructions 
Is  discussed,  and  in  Chapter  IV  passive,  realizable  circuits  are  built  to 
give  the  admittance  of  straight  wire  and  circular  loop  antennas. 


2.4  Sources  of  SEM  Data 


2.4.1  Straight  Wire  Data 


In  a  previous  work  Tesche  [9]  has  derived  the  SEM  data  for  the 
straight  wire,  using  a  method  of  moments  (MoM)  solution  to  Pocklington’s 


Integro-dlfferentlal  equation  for  the  axial  current  flowing  on  the  wire 
I(z,s).  This  equation  has  the  form 


CqE  (z 


.S)  -  f  I(2',8)K(z,z', 


s)  dz'  ,  (2.14) 


where 
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2ir  --Ir 


e  c 


4irR 


ad^ 


(2.15) 


and 


R  -  [(z  -  z’)^  +  d^  sln^  (|)] 


1/2 


(2.16) 


and  L  Is  the  length  of  the  antenna  and  d  Is  the  wire  diameter. 

When  Equation  (2.13)  is  cast  Into  MoM  form,  the  result  is 

^7^  I(ST  -  .  '?.i7) 

Here  z7^  Is  the  system  nxn  matrix,  where  n  Is  the  number  of  zones  on  che 
antenna,  and  I  and  V  are,  respectively,  the  response  and  source  vectors,  each 
of  dimension  n.  Using  this  formulation,  then  the  SEM  poles  are  found  by 
solving 

det  z(s^^)  -  0  (2.18) 


for  the  nontrivial  solutions.  The  natural  current  modes  are  found  from  the 
set  of  equations 

zTsTTT  I(s_^)  -  0  .  (2.19) 
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The  normalization  coefficient  Is  found  by  computing 


^nl 


8  .^nl 
nl 


(2.20) 


Using  these  equations,  an  SEM  data  base  was  constructed  for  a  straight  wire 
with  an  aspect  ratio  (diameter /length)  of  .01.  It  was  found  convenient  to 
use  64  zones  for  the  MoM  equations.  The  resulting  data  agrees  with  the 
results  reported  by  Tesche  [9]. 

The  SEM  poles  for  the  wire  fall  In  layers  In  the  left-hand  s-plane. 

Indicated  In  Figure  d.  Poles  are  indexed  by  (m,n),  where  m  Is  the  layer  and 

n  Is  the  pole,  numbered  sequentially  by  distance  from  the  ~  axis,  where  C  Is 

ire 

the  speed  of  light.  Layer  one  consists  of  those  poles  closest  to  the  ju 
axis,  layer  two  of  the  next  closest,  etc.  Reference  to  this  layer  scheme  of 
ordering  poles  will  be  made  frequently. 

2.4.2  Circular  Loop  Data 

The  SEM  loop  data  Is  an  important  coiq>lement  to  the  numerical  data 
derived  for  the  wire.  In  so  much  as  the  electric  field  Integral  equation  for 
the  loop  can  be  approximately  solved  analytically  due  to  symmetry  by  expand¬ 
ing  the  current  Into  a  Fourier  series,  as  done  by  Wu  [10].  Also  the  loop 
represents  a  doubly-connected  object,  as  opposed  to  the  singly-connected 
wire,  and  the  difference  In  the  SEM  representation  of  the  admittance  for 
these  two  objects  can  be  Investigated.  Namely  the  wire  evidences  zero  admit¬ 
tance  as  s  -*■  0  while  the  loop  admittance  does  not.  Due  to  the  analytical 
tractabllity  of  the  loop  Integral  equation  It  Is  possible  to  Identify  and 
group  the  poles  by  elgenmodes,  which  will  be  useful  when  the  Investigation 
of  elgenadnlttances  Is  made  In  Chapter  III. 
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The  SEM  current  expansion  given  by  Umashankar  [11],  based  on  Wu's 
results.  Is 


I(*,s) 


nir 


z  z 

n  i 


^nl 

®  -  ®ni 


f 


(2.21) 


where  n  Is  the  index  over  modes  and  1  the  index  over  poles  in  a  mode  and  n  is 
the  free  space  Impedance.  V(s)  is  the  transform  of  the  input  waveform.  The 
modal  currents  are  of  the  form  Equation  (2.21)  is  the  representation 

from  which  the  gap  admittance  for  the  loop  is  derived.  Figure  9  gives  the 
loop  geometry. 

When  discussing  the  poles  of  an  eigenmode  for  the  loop,  Wilton  and 
his  co-workers  break  the  poles  into  three  groups  [11,12],  and  this  practice 
will  be  followed  here.  Figure  10  shows  this  classification.  The  type  1  pole 
for  a  mode  is  chat  pole  which  lies  closest  to  the  ju  axis.  Type  II  poles,  of 
which  there  are  a  finite  number,  lie  in  a  semicircular  arc  veering  towards 
the  negative  real  axis.  Type  III  poles,  which  are  infinite  in  number  for  a 
mode,  lie  along  the  Ju  axis.  The  SEM  data  for  the  loop  used  in  this  study 
was  computed  by  Blackburn  [12],  and  has  been  fully  corroborated  by  Umashankar 
[11]. 


CHAPTER  III 


POSITIVE  REAL  CONSIDERATIONS 

3.1  Introduction 

In  this  chapter  the  three  different  modules  (pole-at-a-tlme,  pole- 
pair,  and  terminal  elgenadmlttance)  used  to  construct  the  equivalent  admit¬ 
tance  circuits  are  examined  to  determine  their  PRness  and  hence  their 
realizability.  Admittance  modules  from  both  the  straight  wire  and  circular 
loop  are  examined.  Parameter  tests  for  PRness,  as  defined  In  Chapter  I,  are 
constructed  and  used  extensively.  Also,  the  admittances  of  some  modules  are 
studied  by  graphs  depicting  the  real  part  along  the  Jui  axis.  The  question 
of  the  effect  of  numerical  errors  on  the  PRness  of  modules  Is  addressed. 

3.2  Pole-at-a-Tlme  Circuits 

Ideally,  It  would  be  desirable  If  the  admittance  at  the  feed  point 
could  be  realized  by  constructing  Individual  circuit  modules  on  a  pole  at  a 
time  basis,  as  Illustrated  in  Figure  5.  A  consideration  of  the  Individual 
modules  In  detail  shows  that  such  a  realization  Is  untenable.  The  admittance 
of  an  unmodified  module  In  SEIl  terms  is  given  by 


(3.1) 


This  Is  recognizable  as  an  RL  series  circuit,  with  the  valties  of  the  Inductor 
and  resistor  given  by 

.  _  1 


(3.2) 
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and 

R  -  -  -Si  .  (3.3) 

®nl 

However,  since  the  SEM  data  consists  generally  of  complex  poles  and  residues, 
then  the  L  and  R  elements  take  on  complex  values,  and  are  hence  unrealizable. 
Only  In  the  case  where  a  SEM  pole  lies  on  the  negative  real  axis  Is  the  clr~ 
cult  realizable. 

A  consideration  of  the  modified  pole  at  a  time  circuit  module 
reveals  the  same  situation.  Here  we  have 


(8) 
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(3.4) 


The  general  complex  form  of  the  residue  guarantees  unrealizability. 

Given  that  some  SEM  poles  lie  near  the  ju  axis,  the  question  might 
be  asked:  Will  the  Imaginary  parts  of  the  circuit  elements,  in  either  modi¬ 
fied  or  unmodified  form,  be  negligible?  Taking  the  pole  which  Is  closest  to 
the  jd)  axis  for  the  wire,  for  the  center-driven  case  we  obtain  the  following 
values: 

Unmodified  Modified 

L  -  833.61  -  j233.97  C  --1.36xlo"^  -  jl.32xlo"^ 

R  -  146.39  +  j782.47  R  -  146.39  +  j782.47 

We  can  see  from  this  example  that  the  Imaginary  parts  are  not  negligible. 

3. 3  Conjugate  Pole-Pair  Modules 

From  the  results  In  Section  3.2,  we  see  that  for  realization  of  a 
circuit.  It  Is  necessary  to  have  real  coefficients  for  the  powers  of  s  In  Y  . . 
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This  can  be  done  by  coniiining  conjugate  pole-pairs  Into  one  module.  For  SEM 
data,  this  is  done  as  follows  for  the  unmodified  case 


(s) 
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(3.5) 


where  *  indicates  conjugation.  Simplifying,  we  have 


(8) 


2  Real(a^^)s  -  2  Real(a^s^j^) 
-  2  ReaKs^ps  + 


Similarly  for  the  modified  case  we  have 


(3.6) 


2  Real(a^^s*^)s^  -  2  Real(a*^(s^^)^)s 
|s^^l^(s^  -  2  Real(8j^^)s  +  |Sj^^|^) 


(3.7) 


We  see  that  in  both  cases  the  coefficients  of  s  are  real.  In  addition,  all 
coefficients  in  the  denominators  are  positive,  since  all  poles  lie  in  the 
left-hand  part  of  the  complex  plane. 

In  the  Introduction,  the  conditions  placed  on  the  residues  for  the 
unmodified  admittance  of  Equation  (3.6)  to  be  a  PR  function  are  given.  The 
conditions  for  the  modified  admittance  to  be  PR  are  derived  in  the  following 
paragraphs . 

Rewrite  Equation  (3.7)  as 

d.s^  +  djS 

Y(8)  -  - Y - - -  •  (3.8) 

bj^s  +  b^s  +  b^ 


where 


/ 
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di  -  2  Real(a^^sJ^) 
dj  -  -2  Real(aJ^(s^)^) 

«• 

Then  the  real  part  of  Y(s)  at  s  ■  ju  Is 

d,b-a»^  +  (d,b,  -  d^b  )aj^ 

Real  Y(a)  -  -5—* - - V - 2  *  <3.9) 

b^w^  +  (b2  -  2b^b3)u»^  +  b^ 

Because  the  denominator  formed  In  this  way  Is  the  square  of  an  absolute 
value.  It  Is  always  positive.  Thus  the  question  of  Y(s)  being  positive  real 
devolves  to 

dj^bj^u^  +  (<i2*’2  ■’  ^  u)  i  0  .  (3.10) 

Setting  the  first  derivative  of  this  function  to  zero  to  find  the  extrema 
yields  such  a  point  at  w  *■  0.  Since  this  point  must  be  a  minimum  for  a  PR 
function,  evaluating  the  second  derivative  at  u  -  0  yields  the  first,  low 
frequency,  PR  condition 

2(d2b2  -  d^b^)  >  0  .  (3.11) 

The  other  PR  condition  Is  obtained  from  observing  that  Equation  (3.10),  to  be 

4 

PR  at  high  frequency,  must  have  a  nonnegative  coefficient  for  the  ui  term. 
This  yields  the  condition 

dibi  i  0  . 
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These  conditions  are  seen  to  be  necessary  and  sufficient  for  Equation  (3.7) 
to  be  PR.  Translating  the  coefficients  of  these  two  conditions  Into  the  pole 
and  residue  terns,  where  a^^  -  and  s  ■  we  get 


«nl"nl>  - 


(3.13) 


and 


B  .<*)  . 
nl  nl 


(3.14) 


Equation  (3.13)  Is  called  the  high  frequency  parameter  test  for 
modified  pole-pair  admittances  and  (3.14)  Is  the  low  frequency  parameter 
test,  also  for  modified  pole-pair  admittances.  Both  of  these  tests  must  be 
net  (l.e.,  >  0)  for  the  modified  pole-pair  admittance  to  be  PR.  The  parame¬ 
ter  test  for  unmodified  pole-pair  admittances  Is  derived  in  Chapter  I  and 
repeated  here  as 


"Wni  ■  l®nil“ni  ^  ° 


(3.15) 


For  convenience  we  name  Equation  (3.15)  parameter  test  I,  the  low  frequency 
condition  of  Eqviatlon  (3.14)  parameter  test  IIA,  and  the  high  frequency  con¬ 
dition  of  Equation  (3.13)  parameter  test  IIB. 

The  usefulness  of  these  tests  Is  Illustrated  In  Figures  11  through 
15.  These  figures  display  the  parameter  tests  for  both  modified  and  unmodi¬ 
fied  pole-pair  admittances  for  the  first  five  poles  of  the  first  layer  of  the 
straight  wire,  and  Illustrate  how  the  values  of  the  tests  vary  as  a  function 
of  the  gap  location  on  the  wire.  The  solid  line  Is  the  parameter  test  for 
unmodified  pole-pair  admittances  and  the  dashed  lines  are  for  the  modified 
admittances.  PR  regions  for  both  admittance  forms  are  marked  by  shading. 


Figure  11;  Parameter  tests  for  straight  wire  pole- 
pair  admittances  of  pole  1,  first  layer. 

The  shading  indicates  where  the  PR  conditions 
are  met,  in  either  modified  (dashed  lines) 
or  unmodified  (solid  line)  admittances. 


Figure  13:  Parameter  tests  .for  straight  wire 
pair  admittances  of  pole  3,  first 
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The  x-axls  Is  the  normalized  distance  from  one  end  of  the  antenna,  and  the 
y-axls  Is  the  value  of  the  parameter  tests.  For  the  equivalent  circuit  corre¬ 
sponding  to  a  particular  feed  location  to  be  realizable  on  a  pole-by-pole 
basis,  the  pole-pair  admittance  for  each  pole-pair  must  be  PR. 

Note  that  for  these  first  layer  poles,  which  are  the  primary  con¬ 
tributors  to  the  circuit  admittance  due  to  their  proximity  to  the  ju  axis, 
only  the  laq>ortant  case  of  pole  1  seems  to  be  PR  (In  modified  form)  over  a 
wide  region  of  the  wire.  Also  note  that  In  no  case  for  tlie  first  layer  poles 
is  the  unmodified  form  of  the  admittance  realizable,  and  for  the  modified 
form,  apart  from  pole  1,  only  Isolated  spots  on  the  antenna  appear  PR. 

Figures  16  through  19  Indicate  the  parameter  tests  for  the  first  three  pole- 
pair  admittances  of  the  second  layer  which  lie  off  the  real  axis,  and  the 
second  pole-pair  admittance  of  the  third  layer. 

These  graphs  indicate  only  a  qualitative  measure  of  positive  real¬ 
ness,  however.  If  we  look  at  the  real  part  of  the  modified  conjugate  pole- 
pairs  directly,  greater  Insight  into  the  realizability  Is  available.  Figures 
20  through  27  display  the  real  part  of  the  modified  admittance  along  the  ju 
axis  for  the  first  five  poles  of  the  first  layer  for  the  center  and  quarter 
locations  of  the  wire.  Figures  28  through  32  show  the  real  part  of  the  pole- 
pair  admittances  of  the  second  layer  poles. 

In  these  graphs  the  solid  line  Indicates  the  modified  pole-pair 
admittance  real  part.  The  dashed  line  indicates  the  amount  of  shift  between 
the  modified  and  unmodified  admittance.  In  other  words,  the  dashed  lines 
indicate  the  zero  axis  for  the  modified  form.  Several  facts  can  be  observed 
from  these  graphs.  The  first  layer  poles  are  the  primary  contributors  to  the 
admittance,  as  their  peaks  are  25  to  100  times  greater  than  the  peaks  of  the 


Figure  16;  Parameter  tests  for  straight  wire  pole- 
pair  admittances  of  pole  2,  second  layer 


Figure  17:  Parameter  tests  for  straight  wire  pole- 
pair  admittances  of  pole  3,  second  layer 
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Figure  24:  Real  part  of  straight  wire  pole-pair  ad¬ 
mittance  at  2/L  =  .25  for  pole  2,  first  layer. 
Ther  is  a  slight  negative  area  in  the  low 


Figure  25:  Real  part  of  straight  wire  pole-pair  ad¬ 
mittance  at  z/L  =  .25  for  pole  3,  first  layer 
The  negative  area  is  apparent. 


Figure  26:  Real  part  of  straight  wire  pole-pair  ad¬ 
mittance  at  z/L  =  .25  for  pole  4,  first  layer 
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Figure  28;  Real  part  of  straight  wire  pole-pair  ad¬ 
mittance  at  z/L  =  .5  for  pole  2,  second  layer. 
It  is  entirely  positive. 


Figure  29;  Real  part  of  straight  wire  pole-pair  ad¬ 
mittance  at  z/L  =  .5  for  pole  4,  second  layer 
It  is  entirely  positive. 
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Figure  30:  Real  part  of  straight  wire  pole-pair  ad¬ 
mittance  at  z/L  =  .25  for  pole  2,  second  layer 


32;  Resl  part  of  straight  wir©  pole~pair  ad¬ 
mittance  at  z/L  =  .25  for  pole  4,  second  layer 
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more  deeply  embedded  poles.  The  level  shift  due  to  the  use  of  modified 
admittances  for  first  layer  poles  is  positive,  indicating  that  without  modi¬ 
fication  these  poles  would  all  have  substantial  negative  regions  for  frequen¬ 
cies  below  the  admittance  peaks.  All  of  the  first  layer  poles  except  for 
pole  1  at  the  1/4  point  are  non-PR.  This  result,  as  well  sis  the  PR  quality 
of  pole  1  at  the  1/4  point,  is  predicted  by  the  PR  condition  curves  of  the 
last  section.  However,  now  we  can  see  that  the  violation  of  PRness  is  small 
in  some  cases,  and  consistently  has  the  same  character  for  first  layer  poles; 
l.e.,  the  negative  portion  occurs  as  a  result  of  the  low  frequency  violation 
of  the  PR  conditions  and  results  in  a  small  negative  dip  before  the  main  peak 
of  each  admittance.  This  behavior  suggests  two  questions:  1)  Is  the  small 
negativeness  attributable  to  numerical  error?;  and  2)  If  there  are  no  numer¬ 
ical  errors,  can  the  admittance  still  be  realized  with  perhaps  some  negligible 
negative  components?  It  is  shown  subsequently  that  the  answer  to  question  1 
is  "yes"  only  for  a  few  Isolated  and  predictable  cases.  The  second  of  these 
questions  is  addressed  in  Chapter  IV. 

It  is  worthwhile  at  this  point  to  examine  the  PR  character  of  the 
conjugate  pole-pair  admittances  formed  for  the  circular  loop.  The  type  I 
modified  pole-pairs  for  modes  one  through  five  are  shown  in  Figures  33  through 
37.  These  poles  are  the  major  contributors  to  the  overall  admittance  of  the 
loop.  These  admittances  are  quite  similar  in  character  to  the  first  layer 
pole-pair  admittances  for  the  wire.  Except  for  mode  1,  these  admittances  are 
slightly  non-PR,  exhibiting  a  negative  low  frequency  behavior.  Hie  loop, 
unlike  the  straight  wire,  has  a  nonzero  admittance  at  u  *  0.  This  gives  rise 
to  an  SEM  pole  for  the  loop  at  s  *  0. 


raittance  for  type  I 
is  entirely  positive 


Figure  34:  Real  part  of. circular  loop  pole-pair  ad¬ 
mittance  for  type  I  pole  of  mode  2.  There 
is  a  slight  negative  area  in  the  low  frequency 


Figure  35;  Real  part  of  circular  loop  pole-pair  ad¬ 
mittance  for  type  I  pole  of  mode  3.  There 
is  a  slight  negative  area  in  the  low  frequency 
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Real  part  of  circular  loop  pole-pair  ad¬ 
mittance  for  type  I  pole  of  mode  5.  There 
is  a  slight  negative  area  in  the  low  frequency 


The  character  of  the  modified  pole-pair  admittances  formed  using 
the  type  III  poles  of  mode  2  are  shown  In  Figure  38.  Similar  admittance 
curves  characterize  all  type  III  pole-pair  admittances.  The  characteristics 
are  an  initial  positive  value  rising  to  a  peak,  which  then  falls  to  a  nega¬ 
tive  peak  and  a  negative  final  value  as  frequency  goes  to  Infinity.  The 
peak  value  for  the  first  type  III,  mode  2  pole  Is  240  times  lower  than  the 
peak  of  the  type  I,  mode  2  pole. 

The  modified  pole-pair  admittances  for  the  type  II  poles  of  mode  2 
are  shown  In  Figures  39  and  40.  The  figure  for  the  pole  lying  on  the  nega¬ 
tive  real  axis  Is  not  a  pole-pair,  of  course,  but  only  a  modified  pole  con¬ 
struction.  The  peak  admittance  from  this  pole  Is  some  130  times  lower  than 
the  type  I  pole  of  mode  2.  The  character  of  the  type  II  pole  lying  closest 
to  the  Ju  axis  Is  entirely  negative,  with  a  peak  value  some  50  times  lower 
than  the  type  I  pole  peak.  This  negative  behavior  is  a  characteristic  of  the 
type  II  pole  closest  to  the  ju  axis,  as  Figure  41  indicates. 

In  summary,  the  following  conclusions  regarding  pole-pair  PBness 
are  established: 

-  For  the  straight  wire 

-  Dominant  layer  1  pole-pair  admittances  are  non-PR  in  unmodified 
form,  exhibiting  a  negative  value  at  zero  frequency. 

-  Modified  layer  1  pole-pair  admittances  are  In  general  non-PR, 
exhibiting  only  a  very  slight  negative  value  In  the  low  frequency 
region.  Exceptions  occur  at  isolated  locations  on  the  wire  for 
all  pole-pairs,  and  pole-pair  1  Is  PR  over  most  of  the  antenna. 

-  Deeper  embedded  pole-pairs,  which  contribute  much  less  to  the 
total  admittance,  are  PR  on  various  Intervals,  sometimes  In  modi¬ 
fied  form  and  other  times  unmodified. 

-  For  the  circular  loop 

-  The  dominant  type  I  pole-pairs  are  non-PR  in  unmodified  form, 
exhibiting  a  negative  value  at  zero  frequency. 
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-  Modified  type  1  pole-pairs  are  in  general  non-PR,  e:dilbiting  a 
slight  negative  value  in  the  low  frequency  region.  An  exception 
is  the  pole-pair  of  node  1,  which  is  PR. 

-  Type  II  pole-pairs,  which  contribute  much  less  to  the  total 
admittance  than  do  type  I,  exhibit  a  range  of  PRness.  The  type 
II  pole-pair  closest  to  the  ju  area  is  non-PR  in  unmodified  or 
modified  form,  and  the  pole  admittance  (not  a  pole-pair)  which 
lies  on  the  -a  axis  is  PR.  Other  type  II  pole-pairs  are  not  PR 
in  either  form. 

-  Type  III  pole-pairs,  which  also  have  small  contributions  to  the 
total  admittance,  are  non-PR  in  modified  or  unmodified  form. 

3.4  Eiaenmode  PR  Considerations 

When  the  synthesis  of  equivalent  circuits  is  carried  out  on  a  ter¬ 
minal  eigenadfflittance  basis,  the  theory  of  the  eigenmode  expansion  provides 
some  very  useful  insights  into  the  problem  of  finding  PR  functions.  First  we 
note  that  the  Inverse  eigenvalties  are  positive  real  functions,  as  they 

derive  from  an  impedance  Integral  equation  for  a  passive  object  [6].  The 
PRness  of  the  Inverse  eigenvalues  ,  which  are  termed  elgenadmlttances 

by  Baum  [4],  has  been  demonstrated  by  Wilton  [13].  This  proof  is  reproduced 
in  the  Appendix  for  the  sake  of  completeness.  However,  the  PRness  of  these 
eigenvalues  does  not  Insure  the  PBness  of  the  associated  terminal  elgenadmlt¬ 
tances  of  Equation  (2.13),  as  we  demonstrate  in  the  Appendix.  Additionally, 
the  construction  of  terminal  elgenadmlttances  depends  on  the  availability  of 
the  elgenmode  groupings  of  poles.  For  analytically  tractable  objects,  such 
as  the  loop  or  sphere,  such  groupings  are  available;  however,  for  objects 
such  as  the  wire,  where  the  poles  are  found  numerically,  the  task  of  grouping 
the  poles  by  eigenmodes  is  much  more  difficult.  For  a  complex  object  where 
the  poles  are  experimentally  <.  ived,  if  at  all,  complete  grouping  information 
is  not  likely  to  be  available. 

Wilton  and  his  colleagxies  [13]  have  derived  eigenmode  groups  of 


poles  for  the  wlra  on  a  numerical  basis.  They  have  used  these  groupings  as 
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the  basis  of  conjecture  for  the  general  grouping  scheme  Indicated  In  Figure 
42.  Note  that  there  are  only  a  finite  number  of  poles  belonging  to  each 
eigenmode  In  this  representation.  The  circular  loop  manifests  an  Infinite 
number  of  poles  associated  with  each  eigenmode. 

In  the  case  of  an  Infinite  nunber  of  poles  per  eigenmode,  the 
question  arises  as  to  whether  all  poles  are  necessary  for  PBness  and.  If  not, 
how  many  poles  must  be  Included  to  have  a  PR  function.  In  the  case  of  mode  1, 
where  the  type  I  pole-pair  forms  a  PR  modified  admittance  by  Itself,  Inclu¬ 
sion  of  the  type  II  pole  of  this  mode  destroyed  this  PR  character  In  the  low 
frequency  region,  although  the  maximum  negative  value  Is  only  -.5x10  mhos, 

_3 

compared  to  the  positive  peak  value  of  7.5x10  mhos.  The  systematic  addition 
of  type  III  poles  reduced  the  maximum  negative  value,  and  pushed  the  fre¬ 
quency  at  which  this  maximum  occurred  closer  to  zero  until,  with  the  Inclu¬ 
sion  of  four  type  III  poles,  the  admittance  was  again  PR.  The  Inclusion  of 
these  poles  Is  found  to  have  little  effect  on  the  appearance  of  the  admit¬ 
tance  when  compared  to  the  type  I  pole-pair  admittance,  except  for  introduc¬ 
ing  some  small  high  frequency  variations  due  to  the  type  III  poles.  The  peak 
admittance  value  of  the  eigenmode  grouping  has  been  lowered  by  some  .75  per¬ 
cent. 

An  Interesting  feature  is  observed  when  the  elgenadmlttance  group¬ 
ing  Is  performed.  As  more  type  III  poles  are  included,  the  difference  between 
the  modified  and  unmodified  forms  becomes  smaller,  and  In  the  limit  would 
appear  to  give  the  same  admittance.  In  other  words,  the  constant  terms  Intro¬ 
duced  Into  the  modified  form  sum  to  zero  for  a  given  eigenmode. 


A  similar  examination  of  mode  3  where  the  type  I  pole  Is  not  PR 
revealed  the  same  behavior  as  that  for  mode  1.  Increasing  the  number  of 
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type  III  poles  systematically  reduces  the  negative  portion  of  the  admittance. 

However,  Inclusion  of  all  seven  type  III  poles  available  did  not  produce  a  PR 

—8 

result,  although  the  maximum  negative  value  was  reduced  to  only  -.25>(10  .  For 

higher  modes  a  similar  result  also  holds.  In  general,  to  make  an  elgenadmlt- 
tance  PR  requires  the  Inclusion  of  many  t3rpe  III  poles,  making  the  circuit 
realization  more  complex,  since  the  additional  poles  entail  more  circuit 
elements . 

It  Is  worthwhile  to  point  out  that  If  the  problem  at  hand  were 
merely  to  synthesize  the  total  admittance  of  the  wire  or  loop  using  SEM  data, 
then  the  PR  character  of  the  admittance  Is  assured.  The  mathematical  reason¬ 
ing  leading  to  this  conclusion  Is  as  follows.  If  the  pole  of  the  lowest 
frequency  for  the  structure  is  PR  or  can  be  made  PR,  then  when  all  the  poles 
for  the  structure  are  grouped  the  low  frequency  non-PR  character  of  the 
higher  poles  are  negated  by  the  large  positive  peaks  of  the  lower  poles. 
However,  to  Include  simple  Norton  generators  in  the  synthesis,  it  is  neces¬ 
sary  to  try  to  realize  the  admittance  by  some  group  of  recognizable  modules. 
Further,  It  Is  desirable  that  these  modules  be  as  simple  as  possible. 

When  the  terminal  elgenadmlttances  for  the  straight  wire  are  con¬ 
sidered,  the  results  are  found  to  be  ambiguous.  For  eigenadmlttance  group¬ 
ings  using  conjugate  pole-pairs  at  the  center  of  the  wire,  the  results  are 
PR,  as  Indicated  In  Figures  43  and  44.  However,  when  such  groupings  are 
attempted  for  a  terminal  taken  at  the  quarter  point  of  the  wire,  some  non-PR 
results  are  discovered.  These  groupings  are  Indicated  in  Figures  45  through 
49.  The  postulated  eigenmode  groupings  for  elgenmodes  2,  3,  4,  and  6  are 
oen-PR,  and  In  the  case  of  modes  3  and  6  the  negative  part  Is  not  negligible. 
Tllis  result  could  have  been  foreseen  by  reference  to  the  parameter  test 


Figure  43;  Real  part  of  straight  wire  terminal  eigen- 
admittances  for  mode  3  at  z/L  =  .5.  It  is 
entirely  positive. 


Figure  45:  Real  part  of  straight  wire  terminal  eigen- 
admittance  for  mode  2  at  z/L  =  .25.  There 
is  a  slight  negative  area  in  the  low  frequency 


Figure  47:  Real  p 
admitt 


Figure  49:  Real  part  of  straight  wire  terminal  eigen 
admittance  for  mode  6  at  z/L  =  .25. 
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results.  Taking  the  case  of  node  3,  poles  3  in  the  first  layer  and  2  in  the 
second  layer,  we  see  from  the  FR  condition  graphs  that  at  the  quarter  point 
for  these  poles  the  PR  conditions  for  the  modified  admittance  form  is  vio¬ 
lated  for  both  poles.  Since  the  raal  part  of  each  pole  adds  algebraically 
to  produce  the  total  real  part  admittance,  then  the  grouping  of  non-PR  poles 
cannot  produce  a  PR  result.  Similarly  the  unmodified  admittance  formed  by 
grouping  these  two  poles  is  non-PR. 

While  the  discussion  in  the  Appendix  indicates  that  the  terminal 
elgenadmittances  are  not  demonstratively  PR,  it  is  interesting  to  speculate 
as  to  alternative  reasons  for  the  non-PR  result.  That  numerical  errors  in 
pole/resldue  data  is  the  source  of  departure  from  PBness  is  unlikely  in 
light  of  the  general  error  analysis  discussed  in  the  next  section.  Other 
possibilities  include  Improper  groupings  of  poles  associated  with  the  eigen¬ 
values  and  missing  components  of  the  Mittag-Leffler  expansloti  of  Che  eigen¬ 
values.  Wilton's  groupings  are  based  on  similarity  of  modal  features  and 
the  topological  kinship  of  the  straight  wire  to  the  sphere.  We  are  inclined 
to  trust  Wilton's  conjectured  groupings.  Therefore  it  is  difficult  to  draw 
more  specific  conclusions  at  the  present  time. 

3.5  SEM  Pole/Residue  Error  Effects  on  PR  Considerations 

The  SEM  data  for  the  wire  is  numerically  derived.  This  gives  rise 
to  the  question  as  to  how  strongly  numerical  errors  Influence  the  PRness  of 
the  admittances  derived  from  this  data.  This  is  especially  important  in 
light  of  the  very  small  negative  values  encountered  in  the  modified  conjugate 
pole-pair  admittances  for  first  layer  poles.  We  have  elected  to  investi¬ 
gate  this  question  in  terms  of  the  parameter  tests  given  in  Section  3.2. 
Errors  in  the  real  and  imaginary  parts  of  both  the  poles  and  the  residues 
were  introduced,  and  new  parameter  test  curves  were  produced.  In  such  a  way 
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the  sensitivity  of  these  conditions  to  errors  In  the  poles  and  residues  can 
be  observed.  The  study  is  centered  on  errors  In  the  real  part  of  the  pole. 
The  real  part  of  a  numerically  derived  pole  Is  generally  less  trustworthy 
than  the  Imaginary  part. 

Results  for  pole  1  of  the  wire  are  given  In  Figures  SO  and  51.  In 
these  curves,  the  solid  line  Is  the  original  value,  and  the  dashed  lines 
Indicate  the  changed  values.  Figure  50a  shows  the  effect  of  a  +5  percent 
change  In  the  pole  real  part  on  the  high  frequency  condition  for  pole  1.  The 
change  in  this  condition  Is  very  slight.  Figures  50b  through  50d  Indicate 
the  effect  of  changing  the  other  parameters  by  5  percent.  The  largest  change 
of  some  30  percent,  occurs  when  the  Imaginary  part  of  the  pole  Is  changed. 
Since  this  parameter  is  accurately  known  (within  2  percent)  In  the  numerical 
data,  this  Is  not  bothersome.  Note  that  In  none  of  these  graphs  Is  the  basic 
character  of  the  condition  changed;  chat  is,  the  shape  remains  the  same  and 
the  condition  remains  positive. 

An  entirely  different  situation  occurs  when  we  modify  the  parame¬ 
ters  by  5  percent  for  the  low  frequency  test.  These  effects  are  Illustrated 
In  Figures  51a  and  51b.  In  Figure  51a,  a  change  01+5  percent  in  the  real 
part  of  the  pole  results  In  an  entirely  positive  condition.  This  change 
corresponds  to  forcing  the  real  part  of  the  pole  slightly  away  from  the  ju 
axis.  The  remaining  figures  Indicate  a  similar  occurrence  when  the  Imaginary 
part  of  the  pole  Is  forced  towards  the  real  axis,  when  the  Imaginary  part  of 
the  residue  Is  forced  towards  the  real  axis,  and  when  the  real  part  of  the 
residue  Is  forced  away  from  the  ju  axis. 

The  real  part  of  this  pole  Is  modified  to  produce  a  PR  function, 
both  to  study  the  effect  on  the  admittance  and  for  latter  use  in  circuits. 
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Because  an  error  of  +5  percent  for  the  real  part  of  the  pole  Is  within  the 
computational  uncertainty  of  the  data,  a  value  of  this  pole  with  a  real  part 
changed  by  3.5  percent  Is  used,  which  allows  the  modified  form  of  this  admit¬ 
tance  to  be  PR  over  the  entire  wire.  When  this  Is  done.  It  Is  found  that  the 
peak  value  of  the  modified  admittance  formed  from  this  pole  decreases  by  3 
percent  over  the  original  value. 

In  light  of  these  results  for  pole  1,  the  other  poles  were  studied 
In  the  hope  of  producing  PR  results  for  the  entire  wire  with  adjustments  In 
the  pole  value  within  the  limit  of  numerical  error.  It  was  found  however 
that  all  the  other  poles  retained  the  essential  characteristics  of  the  non-PR 
conditions  under  as  high  an  error  as  10  percent.  Figures  52  through  57  Indi¬ 
cate  these  results  for  poles  2  and  3  for  the  first  layer,  and  pole  2  of  the 
second. 

The  conclusion  to  be  reached  Is  that  the  small  negative  low  fre¬ 
quency  admittance  values  found  when  forming  modified  conjugate  pole-pair 
admittances  are  not  attributable  to  any  small  numerical  errors  In  either  the 
poles  or  residues,  but  are  Inherent  properties  of  these  admittances.  The 
Important  exception  Is  pole  1  In  the  first  layer.  It  was  found  that  when 
large  changes  In  the  poles  or  residues  were  made  In  an  attempt  to  make  these 
pole-pair  admittances  PR,  the  peak  value  of  the  admittance,  which  we  know  to 
be  accurate  through  comparison  with  Integral  equation  admittance  results.  Is 
changed  substantially.  Hence  means  other  than  parameter  modification  or 
approximation  must  be  used  to  realize  these  admittances. 
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(A)  ±10%  CHANGE  OF  POLE  REAL  PART.  (C)  ±107.  CHANGE  OF  RESIDUE  REAL  PART. 

(B)  ±107.  CHANGE  OF  POLE  IMAGINARY  PART.  (0)  ±107.  CHANGE  OF  RESIDUE  IMAGINARY  PART. 

Figure  56:  Effect  of  ±10^  change  on  parameter  test  IIB 
for  pole-pair  admittance  of  pole  2,  second 
layer  of  the  straight  wire. 


CHAPTER  IV 


CIRCUIT  SYNTHESIS 

4.1  Introduction 

This  chapter  addresses  the  construction  of  equivalent  admittance 
circuits  for  three  objects  -  the  center-driven  wire,  the  quarter-driven  wire, 
and  the  circular  loop.  The  circuits  derived  are  physically  realizable;  that 
is,  they  are  composed  only  of  passive  resistor.  Inductor,  and  capacitor 
elements.  The  circuits  are  composed  of  terminal  elgenadmlttance  modules 
whenever  possible.  However,  when  terminal  eigenadmit tances  are  not  PR,  two 
other  methods  are  used  to  create  equivalent  circuits.  In  the  first,  a  method 
is  devised  by  which  an  admittance  having  a  negligible  non-PR  region  can  be 
made  PR.  In  the  second,  poles  are  regrouped  to  produce  a  PR  result.  The 
response  of  these  circuits  to  a  transient  voltage  at  the  port  are  compared  to 
the  response  of  the  original  object  to  the  same  voltage.  A  study  of  the  sen¬ 
sitivity  of  these  circuits  to  element  variations  is  also  made. 

4.2  Center-Driven  Wire 

Two  equivalent  circuits  are  constructed.  The  first  is  composed  of 
terminal  eigenadmlttances,  and  the  second  employs  first  layer  pole-pair  admit¬ 
tances  only.  To  construct  this  second  circuit,  a  method  is  given  by  which 
the  small  non-PRness  of  these  admittances  can  be  neglected.  The  responses  of 
both  circuits  are  investigated. 

4.2.1  Terminal  Elgenadmlttance  Results 

The  PR  results  from  Chapter  III,  which  indicate  that  the  sumnatlon 
of  poles  along  the  elgenmodes  postulated  in  Figure  42  are  PR  and  hence 
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physically  realizable,  are  used  as  a  basis  Co  construct  an  equivalent  circuit 
for  the  center-driven  wire.  The  procedure  followed  Is  to  construct  the 
rational  polynoiiil£LL  formed  by  the  summing  of  those  poles  In  an  eigenmode  and 
their  conjugates,  and  then  use  standard  circuit  synthesis  techniques  to 
derive  the  equivalent  circuits.  The  circuit  synthesis  Is  done  by  reducing 
Che  polynomial  to  a  minimum  conductance  -  minimum  reactance  form,  and  then 
applying  the  Bott-Duffln  technique.  Although  a  complete  representation  for 
the  eigenmode  grouping  Is  available  for  only  the  first  six  eigenmodes,  we 
chose  to  Include  mode  seven,  although  the  data  base  for  this  mode  lacked  one 
pole  In  the  fourth  layer,  since  It  was  PR  without  this  pole.  The  poles  and 
residues  used  are  listed  In  Table  1.  Figure  58  gives  Che  result  of  this 
synthesis,  and  Table  2  lists  the  element  values.  The  reactive  components 
are  normalized  •  The  real  part  of  pole  1  Is  biased  upward  by  3.5  per¬ 

cent  of  Che  peak  value,  as  Indicated  In  Chapter  III,  to  make  It  PR.  Because 
only  modes  which  possess  current  distributions  which  are  even  functions 
couple  at  the  center  of  the  antenna,  eigenmodes  one,  three,  five,  and  seven 
are  realized.  Table  2  also  gives  the  actual  circuit  element  values  needed  to 
realize  the  center-driven  admittance  of  a  100-meter  wire.  This  length  was 
chosen  becatise  It  gives  element  values  which  are  In  the  picofarad,  microhenry 
range.  These  values  change  proportionally  to  length,  so  chat  only  a  certain 
range  of  sizes  for  the  wire  may  actually  be  realizable.  This  scaling  of 
Inductors  and  capacitors  Is  frequency  scaling,  with  the  scaling  factor  equal 
to  ~  .  Hence  both  capacitors  and  inductors  Increase  with  Increasing  length. 


Table  1 

Poles  and  Residues  Used  for  Eigenmode  Synthesis, 
Center-Fed  Cylindrical  Antenna 


Mode 

Poles 

Residues 

1 

-.08427fJ.9158 

.  1112xl0'‘^+j .  3121x10“^ 

3 

-.1473+j2.870 

. 1319x10“ ^+j . 3301x10"^ 

-2.491+J 1.328 

. 0988x10“ . 2529x10"^ 

5 

-,1877+J4.834 

. 1423xl0“^+j . 3521x10“^ 

-2.894+j3.528 

. 1408x10“ .1850x10"^ 

-4.517+jl.497 

. 3218x10“ Vj . 2378x10"^ 

7 

-.2177+J6.792 

. 1496xl0“^+j . 3699x10"^ 

-3.14(>fj5.600 

. 1647x10“ ^+j . 1687x10“^ 

-5.069+j3.890 

. 6818x10“ .1628x10"^ 

Poles  are  normalized  as  per  Tesche  [9] 
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Table  2 

Element  Values  Normalized  to  —  for  Center-Driven  Wire  Admittance, 

ire 

Realized  on  a  Terminal  Elgenadmlttance  Basis 


Mode 


Element  Normalized  Value 


100  Meter  Antenna 
(ohms,  farads,  henries) 


1 

«11 

.663 

^2 

2200.6 

"ll 

2.72xl0‘ 

4l 

434.55 

.663 

2200.6 

288.6x10"^ 

46.11x10“® 


R3I 

29.76 

R22 

4280. 

R33 

571.4 

S4 

53.07 

^35 

397.4 

Si 

290.3 

S2 

3817. 

^33 

141.2 

^34 

18.45 

Si 

3.649x10“^ 

^32 

S3 

8.126x10  ^ 
1.681x10“^ 

S4 

1.278x10  ^ 
-4 

S5 

6.218x10 

29.76 

4280. 

571.4 
53.07 

397.4 


30.80x10  ” 
405.xl0“® 
14.98xl0“® 
1.958x10"® 


38.72x10 


-12 

-12 


8.622x10 
1784x10“^^ 
135.6x10 
65.98x10 


-12 

-12 


Table  2  (continued) 


Element 

Normalized  Value 

100  Meter  Antenna 
(ohms,  farads,  henries) 

*^51 

50.69 

50.69 

R52 

5228 

5228 

.^53 

1003 

1003 

*^54 

2829 

2829 

9A7.2 

9A7.2 

^56 

137.9 

137.9 

^57 

718.5 

718.5 

^58 

761.1 

761.1 

^59 

25A.9 

25A.9 

Si 

1.539x10^ 

16.33x10”^^ 

S2 

3.926x10“^ 

A. 166x10“^^ 

S3 

7.068x10"^ 

7.A99xl0“^^ 

Sa 

3.362x10'*^ 

35.67x10“^^ 

Ss 

2,52xl0~^ 

267.4x10''^^ 

Se 

2.608x10“^ 

27.67x10”^^ 

S7 

1.738x10'"^ 

18.44x10“^^ 

Si 

2A2.A 

25.72x10”® 

S2 

1817. 

192.8x10”® 

S3 

188. 

19.95x10”® 

Sa 

125.3 

13.29x10”® 

Ss 

28.3 

3.003x10“® 

S6 

50.96 

5.407x10”® 

Table  2  (continued) 


Element  Normalized  Value 


100  Meter  Antenna 
(ohms,  farads,  henries) 


^1 

53.02 

^2 

5967 

^73 

2500 

^4 

1682 

^75 

3611 

^6 

106.5 

^7 

254.1 

^78 

377.7 

R79 

175.9 

51 

52 

8.368x10"^ 

2.427x10"® 

53 

54 

Ss 

S6 

2.285x10  ^ 
3.62x10"^ 
2.161x10"® 
4.982x10"^ 

S7 

1.129x10  ^ 

Si 

230. 

S2 

1373 

S3 

316.5 

S4 

717. 

"^75 

15.42 

^76 

14.52 

53.02 

5967 

2500 

1682 

3611 

106.5 

254.1 

377.7 

175.9 

8.879x10"^^ 

2.575x10"^^ 

2.424x10”^^ 

38.41x10'’^^ 

229.3x10“^^ 

52.86x10"^^ 

119.8x10”^^ 

24.40x10"^^ 

145.7x10"® 

33.58x10"® 

76.08x10"® 

1.636x10"® 

1.541x10"® 


i 

J 
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4.2.2  Pole-Pair  Realization  Using  First  Layer  Poles 

Although  the  previous  section  gives  a  practical  circuit  for  the 
admittance  of  a  center-driven  wire,  several  factors  limit  Its  use.  An  eigen¬ 
mode  synthesis  for  modes  higher  than  seven  requires  the  SZM  poles  lying  in  the 
fifth  and  higher  layers.  Such  poles  are  very  hard  to  extract,  and  to  date 
have  not  been  extracted.  This  would  result  in  a  band  limitation  on  the  cir¬ 
cuit  If  the  higher  modes  are  not  represented.  Also,  as  more  modes  are 
Included  In  the  representation  more  poles  per  mode  are  needed  for  realiza¬ 
tion,  which  leads  to  more  circuit  elements  per  mode. 

For  these  reasons  and  because  poles  other  than  first  layer  appear 
to  contribute  negligibly  to  the  total  admittance,  a  circuit  using  only  first 
layer  pole-pair  admittances  Is  constructed.  In  light  of  the  very  small  non- 
PR  values  associated  with  the  first  layer  poles,  it  is  reasonable  to  assume 
that  some  realization  yields  a  circuit  which  Includes  some  small,  negligible 
negative  elements. 

The  following  observation  in  regard  to  the  Bott-Duffin  synthesis 
provided  the  necessary  Insight  for  approximate  realization  of  the  circuit. 
Suppose  a  negative  conductance  G^,  equal  to  the  maximum  negative  value  of  a 
modified  pole-pair  admittance.  Is  removed  from  the  admittance.  The  result¬ 
ing  admittance  Is  PR,  since  its  real  part  has  been  raised  by  a  level  equal  to 
the  absolute  value  of  G^.  If  a  Bott-Duffln  synthesis  Is  then  performed  on 
this  PR  admittance,  the  circuit  module  given  In  Figure  59  results. 

Here  Y  and  Y  are  functions  obtained  in  the  Bott-Duffin  synthesis 

w  Ij 

and  are  siobsequently  synthesized  by  the  removal  of  a  conjugate  pole.  The 
module  Is  seen  to  consist  of  three  branches:  a  shunt  conductance,  a  capaci¬ 
tive  branch,  and  an  Inductive  branch.  Since  the  overall  admittance  of  this 
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module  at  cu  *  0  Is  zero,  then  the  admittance  real  part  of  Y^,  G^,  must  equal 


G^.  Furthermore,  the  conductance  of  the  Inductance  branch  is  limited  to  a 
maximum  of  G  ,  since  the  elements  In  the  branch  are  in  series  with  G  .  Now 

L  la 

since 

iGjll  l^nll 


for  the  first  layer  poles,  then 


I«lI 


«  Y 


cp 

nl 


(5.2) 


as  well. 

This  demonstrates  that  only  the  capacitive  branch  need  be  Included 
in  the  synthesis  and  the  ocher  branches  are  negligible.  This  realization  has 
Che  correct  behavior  at  zero  frequency,  and  has  the  proper  pole  at  the  proper 
peak  value.  The  peak  value  does  not  differ  from  Che  original  because,  at  the 
pole,  Che  inductive  branch  contributes  admittance  G^,  while  the  resistor  con¬ 
tributes  -Gy.  Only  the  final  value  of  the  admittance  would  differ  from  the 
original  pole-pair  admittance,  and  this  is  negligible. 

This  scheme  is  used  Co  synthesize  an  equivalent  circuit  using  only 
first  layer  poles  for  the  center-driven  wire.  The  resulting  circuit  is  given 
in  Figure  60,  and  a  table  of  circuit  element  values  given  in  Table  3. 
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Table  3 


Element  Values  for  Centeir-Drlven  Wire  Equivalent  Circuit 
Realized  by  First  Layer  Pole-Pair  Admittances 


Element 


100  tieter  Antenna 

Normalized  Value  (ohms,  farads,  henries) 


"'ll 

.663 

.663 

^2 

2200.6 

2200.6 

^11 

2.72x10“^ 

288.6x10"^^ 

hi 

434.55 

46.11x10"^ 

hi 

4958.7 

4958.7 

hi 

^22 

4.728x10"^ 

1.063x10”^ 

50.17x10"^^ 

112.8x10"^^ 

hi 

368.22 

39.07x10"^ 

hi 

6513.5 

6513.5 

hi 

C32 

2.253x10"^ 

2.787x10"^ 

-12 

23.9x10 

29.57x10"^^ 

hi 

341.4 

36.22x10"® 

hi 

7542.3 

7542.3 

hi 

h2 

1.418x10“^ 

1.249x10"^ 

15.05x10"^^ 

13.25x10"^^ 

'41 


-6 


324.6 


34.44x10 


Table  3  (continued) 


100  Meter  Antenna 


Pole 

Element 

Normalized  Value 

(ohms,  farads,  henries) 

9 

^51 

8374.8 

8374.8 

Si 

^52 

1.011x10“^ 

7.124x10"^ 

10.73x10"^^ 

7.559x10"^^ 

Si 

311.91 

33.09x10“® 

4.2.3  Circuit  Performance 

The  transient  current  response  at  the  port  of  the  circuits  Is  ana¬ 
lyzed  by  means  of  a  SCEPTRE  clrctilt  analysis  program  Implemented  on  an  IBM 
370  computer.  The  circuit  Is  excited  at  the  port  by  a  Gaussian  pulse  of  the 
form 

V(t)  -  _  ,5  3, 

7  -9 

where  AN  *  3.25^10  and  TMAX  =  60.802x10  seconds.  These  results  are  then 
compared  to  results  from  the  thin  wire  time  domain  (TWTD)  program  for  a  similar 
wire.  Figure  61  Illustrates  these  results  for  the  elgenadmlttance  circuit, 
and  Figure  62  Illustrates  the  results  for  the  first  layer  circuit  synthesis. 

The  response  for  the  two  circuits  Is  almost  Identical.  The  elgen¬ 
admlttance  circuit  shows  more  oscillatory  behavior  but  this  Is  due  to  the 
truncation  of  the  modes  at  mode  seven.  For  late  times  the  response  Is  very 
close  to  the  TWTD  response.  Note  that  both  circuits  miss  the  early  time, 
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Figure  61:  Comparison  of  Che  response  of  Che  eigenadmiCCance 
circuiC  Co  TirTD. 


Comparison  of  Che  response  of  the  first  layer 
pole-pair  circuit  to  TWTD. 
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or  forced,  response.  This  Indicates  chat  the  representation  for  the  wire  Is 
not  complete.  A  study  was  made  to  attempt  to  discover  what  was  missing  to 
cause  this  early  time  response  error.  If  we  observe  that  the  wire  Is  a 
capacitive  object,  then  It  can  be  theorized  that  the  SEM  representation  of 
the  admittance  should  include  a  shunt  capacitor  across  the  port.  Such  a 
capacitor  wotild  correspond  to  a  pole  at  Infinity.  It  can  be  demonstrated 
that  the  response  of  such  a  capacitor  would  be  such  that  It  would  give  the 
correct  forced  response  and  not  affect  the  late  time.  The  current  response 
of  a  capacitor  is  given  by 

l(t)  -  C  ^  .  *  (5.5) 

If  the  forcing  function  Is  a  Gaussian  pulse  then  the  current  contributed  by 
the  shunt  capacitor  Is  given  by 

i(t)  -  C  2(AN(t  -  TMAX))  .  (5.6) 

Such  a  shunt  capacitor  can  be  viewed  as  a  lumping  of  Che  static  capacitances 
associated  with  an  Infinite  number  of  poles  neglected  in  Che  synthesis.  This 
conjecture  is  supported  by  work  done  by  Franceschettl  [14]  on  the  quasl- 
scatlc  capacitance  of  spheroidal  antennas. 

4. 3  Quarter-Driven  Wire 

In  the  last  chapter  we  showed  chat  some  of  the  terminal  elgenadmlt- 
tances  for  this  structure  are  non-PR.  Thus  an  alternate  means  of  realizing  the 
equivalent  circuit  Is  enployed  here.  In  order  to  perform  the  synthesis  at  this 
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feed  point (  we  resorted  to  an  ^  hoc  grouping  of  poles  so  as  to  achieve  PR 
"pole  group"  admittances.  The  groups  used  are  as  follows: 

Group  1 

pole  1,  first  layer 
Group  2 

pole  2,  first  layer 
pole  1,  second  layer 
Group  3 

pole  3,  first  layer 
pole  2,  second  layer 
pole  1,  third  layer 
Group  4 

pole  5,  first  layer 
Group  5 

pole  6,  first  layer 
pole  7,  second  layer 
pole  1,  fourth  layer 
Group  6 

pole  7,  first  layer 
pole  6,  second  layer 
pole  4,  third  layer 
Group  7 

pole  9,  first  layer 
Group  8 

pole  10,  first  layer 
pole  9,  second  layer 
pole  8,  third  layer 


-  104  - 


Poles  4  and  8  of  the  first  layer  were  not  Included,  because  of  their  negligi¬ 
ble  pole-pair  admittances.  Table  4  lists  the  values  of  the  poles  and  residues 
used. 


Table  4 

Poles  and  Residues  Used  for  1/4-Fed  Synthesis  of  Cylindrical  Antenna 


ad  hoc*  Group 
1 
2 

3 

4 

5 

6 

7 

8 


Poles 

-. 0842 7+j. 9158 

-.1199+jl.890 

-2.149+jO.O 

-.1473^-j2.870 

-2.491+jl.328 

-4.098+J0.0 

-.1877+j4.834 

-.2038fj5.814 
-3.2254-j6.620 
-6. 006- j.  884x10' 

-.21774-j6.792 

-3.297-t-J7.632 

-5.069+j3.890 

-.24264-j8.736 

-.1783+j9.766 

-3.363fj8.636 

-5.6384-j8.235 


Poles  are  normalized  as  per  Tesche  [9] 
*Resldue  value  unavailable 


Residues 


. 6192x10“ ^+j. 1619x10' 

.1268xl0“^4-j.  3537x10' 
-.  1707x10"  ^-t-j.O 

.5496xl0"^+j. 2259x10' 
.1159xl0"^-j.  1483x10' 
-.2497xl0“\j0.0 

.  9177x10"  ^+j.  1117x10' 

.  1500x10"  ^4-j.  4254x10' 
.  9804x10"  ^-j .  2504x10' 

*  * 


.5494xl0"^+j . 3144x10" 
. 5421x10" Vj . 1116x10"' 
.  6790x10"  ^-l-j .  8824x10" 

.  1074xlo"^4-j .  7253x10"' 

.  1641xl0"^4-j .  5373x10" 
-.  1079x10"  ^-t-j .  6414x10' 
.2912xl0"^+j.8177xl0" 


Figure  63  gives  the  circuit  synthesized  from  the  above  groupings 
using  the  Bott-Duffln  procedure.  Table  5  gives  the  element  values.  The  cir¬ 
cuit  was  analyzed  with  the  SCEPTRE  program,  using  the  same  Gaussian  pulse  as 
In  the  previous  example,  and  a  con^arlson  with  results  computed  by  TWTD  was 
made.  Figure  6A  Illustrates  this  comparison.  The  circuit  response  exhibited 
(when  compared  to  TWTD)  oscillations  of  relatively  high  magnitude.  These 
oscillations  cause  the  signal  to  be  degraded  In  late  time.  However,  the  gen¬ 
eral  shape  was  replicated.  The  early  time  response  shows  the  same  departure 
from  TWTD  as  did  the  center-driven  wire.  There  are  two  possible  sources  of 
the  oscillations.  One  Is  an  Incorrect  element  value  in  one  of  the  modules. 

The  other  more  likely  possibility  is  an  error  in  the  SEM  poles,  partic¬ 
ularly  In  the  higher  frequency,  deeply  endiedded  ones. 

4.4  Admittance  Synthesis  for  Circular  Loop 

As  stated  In  Chapter  III,  there  are  an  infinite  number  of  poles  in 
an  eigenmode  of  the  circular  loop*  Therefore,  a  terminal  eigenadmlttance 
module  has  an  Infinite  number  of  elements  and  Is  unrealizable  in  a  practical 
sense.  By  truncating  the  set  of  eigenmode  poles,  a  truncated  terminal  elgen- 
admlttance  module  can  be  constructed.  Two  problems  remain,  however.  First, 
the  tmncatlon  may  not  result  In  a  PR  admittance.  In  which  case  the  modified 
Bott-Duffln  procedure  of  Section  4.2.2  must  be  used.  Namely,  the  admittance 
Is  made  PR  by  removing  a  negative  conductance,  and  then  neglecting  this  con¬ 
ductance  branch  and  the  Inductance  branch  In  the  subsequent  Bott-Duffin 
synthesis.  Second,  even  If  the  truncated  admittance  Is  PR,  the  complexity  of 

the  Bott-Duffln  circuit  grows  quickly  with  pole  count.  The  number  of  elements 

N 

required  for  a  Bott-Duffln  synthesis  grows  by  7x2  -6,  where  N  Is  the  number  of 
pole-pairs.  Thus  an  admittance  composed  of  4  pole-pairs  requires  106  elements. 


63:  Equivalent  circuit  for  quarter-driven  antenna 


/ 
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Table  5 


Pole 

Element  Values 

for  Pole  Groupings  of 

Quarter-Driven  Antenna 

100  Meter  Antenna 

Group 

Element 

Normalized  Value 

(ohms,  farads,  henries) 

1 

\l 

7.427 

7.427 

^2 

4311.9 

4311.9 

"ll 

1.51x10“^ 

160.2x10“^^ 

-6 

hi 

782.3 

83x10 

2 

hi 

2234.4 

2234.4 

h2 

3.405 

3.405 

ha 

.752 

.752 

C-, 

a.ooaxio"*^ 

84.91x10“^^ 

21 

c„- 

2.288x10'*^ 

24.28xl0‘^^ 

22 

1.787x10“^ 

18.96x10“^ 

23 

h4 

1.026x10^ 

108.9x10"^ 

-6 

hi 

300.5 

31.88x10 

-3 

^22 

17256. 

1.831x10 

-9 

ha 

.3846 

40.81x10  ^ 

3 

hi 

36.13 

36.13 

R32 

5806. 

5806. 

ha 

4677 

4677 

h4 

1597 

1597 

hs 

1042 

1042 

he 

355.6 

355.6 

h7 

286.5 

286.5 

/ 
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Table  5  (continued) 


Pole 

Group 

Element 

Normalized  Value 

100  Meter  Antenna 
(ohms,  farads,  henries) 

3 

Si 

1.843x10"^ 

-5 

19.56x10'"^ 

-12 

C,o 

6.606x10 

7.009x10 

32 

-5 

-12 

3.24x10  ^ 

3.438x10 

33 

C-,, 

3.174x10"^ 

33.68x10“"^ 

34 

1.634x10"^ 

17.34x10“"^ 

35 

-3 

-12 

Se 

2.007x10 

213.0x10 

527.9 

56.02x10”^ 

31 

^32 

3339 

354.2x10"*’ 

-6 

271.8 

28.84x10  ” 

33 

"34 

107.9 

11.66x10"^ 

--6 

"35 

53.9 

5.719x10 

4 

Does  not 

couple 

5 

Si 

1.418 

1.418 

R52 

32510. 

32510. 

Si 

7.893x10"** 

837.5x10"^^ 

"51 

541.3 

57.43x10"^ 

6 

Si 

65.91 

65.91 

S2 

17640. 

17640. 

S3 

38400. 

38400. 

S4 

32450. 

32450. 

S5 

258.3 

258.3 

Se 

118.7 

118.7 

S7 

140.4 

140.4 
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Table  5  (continued) 


Pole 

Group 


100  Meter  Antenna 

Element  Normalized  Value  (ohms,  farads,  henries) 


Si 

1.182x10"^ 

12.5Axl0"^^ 

Sa 

i.assxio"^ 

1.335x10"^^ 

S3 

1. 790x10”® 

1899x10"^^ 

Sa 

A. 909x10"® 

5.209x10"^^ 

Ss 

7.12Axl0"^ 

75.58x10"^^ 

Se 

2.A71xl0"^ 

26.22x10"^^ 

Si 

223.7 

23.7Axl0"® 

Sa 

32  A6. 

3AA.Axl0"® 

S3 

1126. 

119.5x10"® 

Sa 

57.33 

6.083x10"® 

Ss 

8.156 

.865Axl0"® 

Si 

108.3 

108.3 

Sa 

2615. 

2615. 

S3 

708. 

708. 

Sa 

65.A 

65. A 

Ss 

31.12 

31.12 

Se 

336.9 

336.9 

S7 

2661. 

2661. 

00 

28810. 

28810. 

S9 

720.5 

720.5 

SlO 

605A0 

605A0 

^711 

5592. 

5592. 
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Pole 

Group 

7 


Table  5  (continued) 


Element 


Normalized  Value 


100  Meter  Antenna 
(ohms,  farads,  henries) 


'71 

72 

73 

74 

75 

76 

77 

78 

79 

'710 

'711 

'712 


9.245x10 

1.065x10"^ 

2.084x10"^ 

3.532x10“^ 

1.642x10"^ 

3.355x10"^ 

8.445x10"^ 

1.277x10"^ 

4.131x10"^ 

1.921x10'"^ 

3.924x10"® 

2.437x10"® 


L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 


71 

72 

73 

74 

75 

76 

77 

78 

79 

710 

711 


159.1 

240.6 
77.82 
361.9 
7.392 
4.592 

200.6 
392.6 
6654. 
30940. 
632.0 


-12 


9.809x10 

.-12 


11.3x10 

22.11x10 


-12 


374.8x10 

,-12 


1743x10 

35.6x10' 


-12 


8.961x10 

13.55x10 

4.383x10 

20.38x10 

.4163x10 


-12 

-12 

ri2 

-12 

-12 


.2586x10 


,-12 


16.88x10 

25.52x10"® 

8.257x10"® 

38.39x10"® 

. 7843x10"® 

.4872x10"® 

21.28x10"® 

41.66x10"® 

706.0x10"® 

3.283x10"^ 

67.06x10"® 


8 


Does  not  couple 


-  Ill  - 


Table  S  (continued) 


Pole 

Group 


10 


1. 

i 


Rlement 


Normalized  Value 


^1 

64.38 

R92 

102300. 

Si 

2.819x10"' 

Si 

464.2 

S.01 

5.875 

S02 

3920. 

So  3 

4209. 

S04 

1610.0 

S05 

694.0 

S06 

•  646.3 

S07 

84.55 

*^108 

80150. 

S09 

32170. 

SolO 

1690. 

Soil 

33.94 

^101 

6.284x10"' 

^102 

2.144x10" 

S03 

5.458x10" 

1 

^104 

4.066x10 

S05 

2.147x10" 

^106 

2.418x10" 

So  7 

5.348x10" 

^108 

1.013x10" 

^109 

1.359x10" 

^1010 

4.989x10" 

Soil 

7. 346x10" 

S012 

1.830x10" 

100  Meter  Antenna 
(olins,  farads,  henries) 

64.38 

102300. 


2.991’‘10 

49.26x10 

5.875 
3920. 
4209. 
1610. 
694. 
646.3 
84.55 
80150. 

12 1  70 
1690. 
33.94 

6.668x10 

2.275x10 

5.791x10 

4.314x10 

2.278x10 

25.66x10 

5.674x10 

10.75x10 

14.42x10 

5.294x10 

77.94x10 

19.42x10 


-12 


-6 


-12 

-12 

-12 

-12 

,-12 

-12 

-12 

1 

-12 

I 

-12 

-12 

-12 
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Pole 

Group 

10 


Table  5  (continued) 


100  Meter  Antenna 

Element  Normalized  Value  (ohms,  farads,  henries) 


‘'lOl 

402 

403 

^104 

^05 

^106 

^107 

^108 

^109 

"1010 

"1011 


135.7 
658. 

275.5 

145.5 

369.8 
58.33 

148.4 

58.4 

110.7 

497.8 
1998. 


14.40x10  ” 
69.82x10“^ 
29.23x10"® 
15.44x10"® 
39.24x10"® 
6.189x10"® 
15.75x10"® 
6.196x10"® 
11. 74x10"® 
52.82x10"® 
212.0x10"® 
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For  these  reasons »  two  circuits  for  the  circular  loop  are  con¬ 
structed,  representing  two  degrees  of  complexity.  In  the  first,  only  the 
type  1  poles  of  each  elgenmode  are  used.  In  the  second,  all  poles  necessary 
for  a  PR  result  are  used,  except  when  the  data  Is  Insufficient  for  a  PR 
result.  In  which  case  the  type  I  and  type  III  poles  are  xjsed. 

4,4.1  Type  I  Poles 

Given  that  the  tjrpe  I  poles  are  the  major  contributors  to  the  total 
admittance  of  the  loop,  and  that  the  non-PR  excursions  of  these  pole-pair 
admittances  are  quite  small,  a  modified  Bott-Duffln  synthesis  of  the  type 
used  for  realizing  the  first  layer  poles  on  the  straight  wire  Is  Implemented. 
Because  the  loop  Is  a  closed  object.  It  has  a  pole  located  at  s  >  0.  This 
pole  represents  the  magnetostatic  Inductance  of  the  loop  which  was  realized 
straightforwardly  as  a  lumped  Inductor.  Only  the  type  I  poles  for  the  first 
10  modes  were  used  In  the  equivalent  circuit.  These  are  sufficient  to  real¬ 
ize  the  response  to  the  bandwidth  of  excitation  used.  The  resulting  circuit 
Is  shown  in  Figure  65.  Table  6  lists  the  element  values,  both  normalized  and 
for  a  100  meter  radius  antenna. 

Again  the  circuit  was  analyzed  by  SCEPTRE  and  compared  to  TWTD 
results.  Figure  66  displays  this  comparison.  The  agreement  in  early  time  is 
excellent,  but  the  peak  values  of  the  circuit  are  approximately  15  percent 
higher  than  the  TWTD  response.  In  the  late  time  the  circuit  response  Is  off¬ 
set  from  the  TWTD  response  by  about  0.2  ma,  and  the  circuit  response  shows 


oscillations 
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Table  6 

Element  Values  for  Circular  Loop  Equivalent  Circuit  Using  Type  I  Poles 


Mode 

Element 

Normalized  Value 

100  Meter  Antenna 
(ohms,  farads,  henries) 

0 

^01 

2165. 

721.5x10"^ 

1 

^11 

47.35 

47.35 

«12 

10.16x10^ 

10.16x10^ 

"ll 

9.918x10"^ 

330.6x10“^^ 

4l 

925.5 

308.5x10“® 

2 

^21 

16.63x10^ 

16.63x10^ 

Si 

2.802x10“^ 

93.4x10“^^ 

S2 

1.335x10"^ 

4450. xio“^^ 

Si 

862.5 

287. 5> 

3 

Si 

20.71x10^ 

20.71x10^ 

Si 

1.529x10"^ 

50.98x10“^^ 

S2 

8.616x10"^ 

287.1x10“^^ 

Si 

818.1 

272.7x10“® 

4 

Si 

24.21x10^ 

24.21x10^ 

Si 

1.007x10"^ 

33.58x10“^^ 

S2 

3.189x10"^ 

106.3x10“^^ 

Si 

785.7 

261.9x10“® 
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Table  6  (continued) 


Element 

Normalized  Value 

100  Meter  Antenna 
(ohms,  farads,  henries) 

^51 

26.98x10^ 

26.98x10^ 

51 

52 

7.332x10“^ 

1.660x10“^ 

24.44x10"^^ 

55.32x10“^^ 

Si 

760.2 

253.4x10"^ 

Si 

29.24x10^ 

29.24x10^ 

51 

52 

5.682x10^ 

1.016x10"^ 

18.94x10"^^ 

33.86x10"^^ 

Si 

738.9 

246.3x10"^ 

Si 

31.11x10^ 

31.11x10^ 

Si 

^72 

4.596x10"^ 

6.858x10"^ 

15.32x10"^^ 

22.86x10"^^ 

Si 

720.9 

240.3x10"® 

Si 

32.68x10^ 

32.68x10^ 

Si 

^82 

3.831xl0'^ 

4.938x10“^ 

12.77x10”^^ 

16.46x10"^^ 

Si 

705.3 

235.1x10"® 

I 
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Table  6  (continued) 


Element 

Normalized  Value 

100  Meter  Antenna 
(ohms,  farads,  henries) 

^1 

34.02x10^ 

34.02x10^ 

Si 

3.27x10“^ 

10.9x10“^^ 

S2 

3.729x10"^ 

12.43x10"^^ 

Si 

691.2 

230.4x10“® 

Soi 

35.17x10^ 

35.17x10^ 

Soi 

2.838x10“^ 

9.46x10“^^ 

^102 

2.916x10“^ 

9.72x10“^^ 

226.2x10 


678.6 


Figure  66:  Comparison  of  the  response  of  the  type 
pole-pair  circuit  to  TWTD. 


4.4.2  Type  I  and  Type  III  Poles 


The  departure  of  the  transient  current  response  of  the  circuit 
from  that  of  the  TWTD  model  in  late  time  displayed  in  Figure  65  for  the  type 
1  pole-pair  synthesis  is  attributable  to  the  omission  of  type  III  poles  in 
the  frequency  range  of  Interest.  Therefore  a  circuit  is  constructed  vdilch 
Includes  all  type  III  poles  for  each  eigenmode,  up  to  the  cutoff  frequency 
given  by  mode  10.  For  mode  0,  the  type  II  pole  is  grouped  vlth  the  type  III 
poles,  since  this  results  in  a  PR  function.  For  modes  1  and  2  and  type  I 
plus  type  III  pole  groupings,  the  admittances  are  PR;  for  higher  modes,  they 
are  not.  For  these  higher  modes  the  modified  Bott-Duffln  procedure  is  used 
to  make  the  groupings  PR.  The  circuit  derived  is  given  in  Figure  67,  with 
element  values  given  in  Table  7. 

An  analysis  of  this  circuit  on  SCEPTRE  and  comparison  to  TWTD 
showed  an  improved  response,  particularly  in  the  late  time,  as  indicated  in 
Figure  68. 

4.5  Sensitivity  of  Circuits 

Although  the  circuits  derived  for  the  straight  wire  and  the  loop 
are  in  principle  realizable,  practicality  of  the  realization  still  may  be 
limited  by  sensitivity  considerations.  In  order  to  assess  the  sensitivity, 
the  two  circuits  previously  described  for  the  center-driven  wire  are  sub¬ 
jected  to  pseudorandom  changes  in  element  value  over  a  range  of  +10%  to  -10%. 
The  random  number  generator  la  given  in  Reference  [19].  In  the  range  of  -10 
to  +10  it  has  a  mean  of  zero  and  a  standard  deviation  of  5.8.  No  attempt 
was  made  to  do  a  complete  Monte  Carlo  or  worst- case  analysis  on  circuit  per¬ 
formance.  Several  circuits  with  errors  were  made.  The  resulting  circuits 


Lgure  67a:  E  lulvalent  circuit  for  the  circular  loop  using 
type  I  and  type  III  poles,  modes  0  and  1. 


/ 


Mode 

0 
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Table  7 

Element  Values  for  Equivalent  Circuit  for  the  Circular 
Loop  Using  Type  I  and  III  Poles 


100  Meter  Antenna 


Element 

Normalized  Value 

(ohms,  farads,  henries) 

«01 

12.6x10^ 

12.6x10^ 

^02 

22.34x10^ 

22.34x10^ 

^03 

1.123x10^ 

1.123x10^ 

^04 

120.7 

120.7 

^05 

819.2 

819.2 

^06 

408.4x10^ 

408.4x10^ 

^)7 

1143 

1143. 

®08 

122.9 

122.9 

^09 

168.4 

168.4 

^010 

401.3x10^ 

401. 3x10^ 

®D11 

2208 

A 

2208 

®012 

60.2x10^ 

60.2x10^ 

®013 

43.92x10^ 

43.92x10^ 

^014 

43.15x10^ 

43.15x10^ 

^015 

292.8x10^ 

292.8x10^ 

"oi 

5.334x10"^ 

17.78x10“^^ 

S2 

1.233x10"^ 

4.111x10"^^ 

S3 

8.289x10"^ 

2.763x10"^^ 

V 

2.473x10"^ 

8.242x10"^^ 

>=05 

2.896x10"^ 

965.3x10"^^ 

Se 

1.089x10“^ 

3.631x10"^^ 

S7 

1.699x10"^ 

5663x10"^^ 

Ss 

1.520xl0"^ 

5.065x10"^^ 

S9 

6.999x10"^ 

23330x10"^^ 

SlO 

2.076x10"^ 

691.9x10"^^ 

Sii 

4.752x10"^ 

15.84x10"^^ 

Table  7  (continued) 


Element 

Normalized  Value 

100  Meter  Antenna 
(ohms,  farads,  henries) 

Sl2 

1.772x10“^ 

5.907x10"^^ 

Sl3 

2.319x10"® 

.00773x10"^^ 

Sl4 

6.486x10"^ 

21.62x10"^^ 

^015 

3.048x10"® 

.01016x10"^^ 

Sl6 

1.958x10"^ 

65.27x10"^^ 

Sl7 

6.915x10"® 

.02305x10"^^ 

Sl8 

4.260x10"® 

.0142x10"^^ 

^019 

1.552x10"^ 

51.73x10"^^ 

S20 

5.808x10"® 

1.936x10"^^ 

S2I 

4.959x10"® 

.01653x10"^^ 

^022 

8.103x10"® 

2.701x10"^^ 

^00 

2165 

721.5x10"® 

^01 

7653 

2551x10"® 

^02 

9657 

3219x10"® 

^03 

286.4 

95.46x10"® 

^04 

2344 

781.3x10"® 

^05 

1.503 

.5009x10"® 

"06 

2.445 

.815x10"® 

^07 

3198 

1066x10"® 

^08 

1.144 

.3812x10"® 

^09 

3.411 

1.132x10"® 

SlO 

399.6 

133.2x10“® 

"on 

2.1 

. 7x10"® 

"012 

8.376x10® 

.1792 

"013 

408.6 

136.2x10"® 

"014 

608.1 

202.7x10"® 

"015 

1219. 

406.4x10"® 

"016 

1.428x10® 

47598x10"® 

Table  7  (continued) 


100  Meter  Antenna 

Element  Normalized  Value  (ohms,  farads,  henries) 


"017 

^018 

"019 

"020 

"021 


\l 

^2 

^13 

«15 

^7 

^8 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


11 

12 

13 

14 

15 

16 

17 

18 
19 
110 


1.026x10^ 

.03412 

749.4 

249.8x10"® 

537.0 

179. xio"® 

873.9 

291.3x10"® 

3.45x10® 

1.15 

49.36 

49.36 

1182 

1182 

27.07x10^ 

27.07x10® 

2294 

2294 

8.314 

8.314 

111.6x10^ 

111.6x10® 

9460 

9460 

30801x10^ 

30801x10® 

9.891x10"^ 

329.7x10"®^ 

6.309x10"® 

2.103x10"^^ 

1.958x10"® 

6.528x10"^^ 

8.027x10"^ 

26757. xio"^® 

1.258 

.4193x10"® 

5.979x10"® 

1.993x10"^® 

6.636x10"® 

2.212x10"^^ 

9.369x10"® 

31.23x10"^® 

8.433x10"® 

2.811x10"^^ 

1.459x10"^ 

.4862x10'^® 

126  - 


Table  7  (continued) 


Element 


Normalized  Value 


100  Meter  Antenna 
(ohms,  farads,  henries) 


hi 

1700 

566.5x10"® 

h2 

2.399x10^ 

7997x10"® 

h3 

1531 

510. 3x10"® 

h4 

2159 

719.7x10"® 

hs 

930.3 

310.1x10"® 

h6 

.3735 

.1245x10"® 

h7 

20.56x10® 

6.852 

hs 

1616 

538.6x10"® 

h9 

5016 

1672x10"® 

^110 

322.2x10® 

107.4 

hi 

1.521 

1.521 

h2 

5199. 

5199. 

ha 

29.63x10^ 

29.63x10^ 

h4 

12.69x10^ 

12.69x10^ 

CM 

2.734x10“^ 

91.14x10"^^ 

h2 

6.549x10"® 

2.183x10"^^ 

ha 

6.627x10"® 

2.209x10"^^ 

h4 

7.314x10"® 

2.438x10"^^ 

hi 

2750 

916.5x10"® 

^22 

2492 

830.6x10"® 

ha 

868.5 

289.5x10"® 

CM 

2462 

820.8x10"® 

127 


Element 


Table  7  (continued) 


Normalized  Value 


62.01x10- 


181.9x10- 

27.98x10- 

1.469x10' 

2.535x10“ 

1.090x10" 

2.242x10" 

2.448x10" 


823.5 

4398. 


71.58x10- 

212.4x10' 

34.31x10' 


9,588x10 

2.046x10" 

3.735x10" 

1.607x10" 

1.644x10" 


789.9 


100  Meter  Antenna 
(ohms,  farads,  henries) 

62.01x10^ 

181.9x10^ 

27.98x10^ 


48.95x10 

,8451x10" 

363.4x10" 

.7474x10" 

.8159x10" 


1416x10  ° 
1297. xio"^ 
274.5x10"^ 
1466x10”® 


71.58x10- 

212.4x10' 

34.31x10' 


31.96x10  ' 
.6821x10"^' 
124.5x10"^' 
.5355x10"^" 
.548x10"^^ 


1346x10  ' 
1315x10"* 
263.3x10" 
1675x10"* 


-J. 


-  128  - 


Table  7  (continued) 


Mode 

5 


Element 


Normalized  Value 


iiecer  Antenna 
(ohms,  farads,  henries) 


*^51 

79.92x10^ 

R 

79.92x10^ 

52 

215.6xio'^ 

«53 

40. 3x10^ 

215.6x10^ 

40.3x10^ 

Si 

6.942x10~^ 

23.14x10“^^ 

c 

S2 

1.732x10  ® 

.5773x10“^^ 

62.78xl0“^2 

.4095x10“^^ 

S3 

1.883x10“^ 

S4 

c 

1.229x10“® 

Ss 

1.151x10  ® 

.3837x10“^^ 

Si 

3804 

L 

3957 

1268x10  ° 

52 

L 

763.2 

1319x10  ® 

53 

L 

5580 

254.4x10  ® 

S4 

1860x10“® 

Si 

29.24x10^ 

29.24x10^ 

Si 

5.682x10“^ 

18.94x10“^^ 

c 

62 

1.016x10  ^ 

33.86x10“^^ 

Si 

738.9 

246. 3xlo“® 

Si 

31.11x10^ 

31.11x10^ 

Si 

4.596x10“^ 

15.32x1o“^2 

c 

ni 

6.858x10  ^ 

22.86x1o“^2 

'71 


240.3x10“® 


720.9 
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/ 


Table  7  (continued) 


100  Meter  Antenna 


Element 

Normalized  Value 

(ohms,  farads,  henries) 

32.68x10^ 

32.68x10^ 

^81 

^82 

3.831x10“^ 

4.938x10“^ 

12.77x10"^^ 

16.46x10"^^ 

00 

705.3 

235.1x10"^ 

Si 

34.02x10^ 

34.02x10^ 

51 

52 

3.27x10"^ 

3.729x10"^ 

10.9x10“^^ 

12. 43x10*’^^ 

Si 

691.2 

230.4x10"^ 

Soi 

35.17x10^ 

35.17x10^ 

^101 

^102 

2.838x10"^ 

2.916x10"^ 

9.46x10"^^ 

9.72x10"^^ 

^101 

678.6 

226.2xlo”® 
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were  then  run  on  SCEPTRE.  A  typical  result,  compared  to  the  uncorrupted  run, 
is  given  In  Figure  69  for  the  first  layer  pole-pair  synthesis.  The  responses 
are  Identical  In  the  early  time,  where  only  the  object's  force  response  Is 
Important.  In  late  time,  the  pole  shifting  caused  by  the  element  variations 
results  In  a  distorted  response.  The  sensitivity  of  the  elgenadmlttance 
circuit  was  also  examined  this  way,  with  a  similar  result.  This  Is  not  sur¬ 
prising,  since  It  Is  well  known  that  a  Bott-Duffln  synthesis  Is  very  sensitive 
to  element  errors  [15]. 


w 


CHAPTER  V 


CCRICLUSIONS 


We  have  denonstrated  the  feasibility  of  synthesizing  passive^ 
component  equivalent  circuits  which  simulate  a  single-port  innnlttance  of  a 
passive  antenna/ scatterer.  Circuits  developed  here  replicate  well  the 
terminal  response  of  straight-wire  and  wire-loop  structures  except  In  the 
early  time.  It  Is  conjectured  that  clrciilts  which  account  for  the  static 
capacitance  of  the  structure  completely^  can  bring  this  early  time  behavior 
Into  agreement,  as  well. 

From  the  point  of  view  of  developing  simple  systematic  equivalent 
circuits,  one  would  desire  that  pole-pair  contributions  to  Immlttance  would 
prove  to  be  positive  real  functions.  Disappointingly,  this  study  reveals  an 
abundance  of  counterexamples  to  this  desirable  result.  We  have  observed, 
however,  that  the  dominant  poles  In  the  resonant  structure  of  the  object  - 
those  nearest  the  ju  axis  -  have  admittances  which  manifest  a  near  PR  behav¬ 
ior.  Through  negligible  adjustments  to  these  admittances,  PR  behavior  was 
achieved. 

Terminal  eigenadmlttances  for  the  circular  loop  and  center-driven 
wire  are  shown  to  be  PR  within  the  limits  of  numerical  accuracy.  However, 
certain  terminal  eigenadmlttances  for  the  quarter-driven  wire  are  non-PR. 
Althou^  It  Is  known  that  the  Inverse  eigenvalues  are  PR,  the  PRness  of  all 
terminal  eigenadmlttances  remains  to  be  demonstrated.  Since  the  elgenset 
grouping  of  poles  Is  only  conjectured.  It  Is  not  possible  to  reach  a  definite 
conclusion  at  this  time  In  regard  to  terminal  elgenadmlttance  PRness. 
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The  spherical  antenna  Is  suggested  as  an  object  for  study  of  admit¬ 
tance  synthesis  properties  In  the  future.  That  elgensets  are  well-defined 
and  that  each  elgenset  embraces  a  finite  collection  of  poles  provide  the 
potential  for  complete  and  conclusive  studies. 

The  present  work  has  addressed  only  the  admittance  element  of 
Norton  equivalent  modules.  The  companion  source  terms  need  to  be  studied 
In  light  of  the  present  work.  In  particular,  when  sources  are  combined  even 
on  a  pole-pair  basis,  the  composite  source  exhibits  a  frequency  dependent 
transfer  function  from  the  excitation  waveform.  This  frequency  dependence 
Introduces  circuit  complexities  which  are  likely  to  preclude  physical  real¬ 
ization.  Approximate  methods  for  circumventing  this  frequency  dependence 


are  warranted 


/ 


APPENDIX 


The  proof  of  the  positive  realness  of  the  Integral  operators  for 
the  antenna  scattering  problem  Is  given  here.  This  proof  Is  due  to  Wilton 
[13]  and  reproduced  here. 

Consider  a  conducting  scatterer  S  with  unit  surface  normal  n  lllu^ 

mlnated  by  an  Incident  field  A  surface  current  J  Is  Induced  on  S 

“S  — S 

fdilch  produces  a  scattered  field  (E  ,H  )  satisfying  the  boundary  condition 

n  X  (E^  +  E^)  =  0  (A.l) 

on  S.  By  the  equivalence  theorem,  the  scatterer  may  be  removed  and  replaced 
by  the  surface  current  J  radiating  In  free  space.  This  current,  radiating 
In  the  absence  of  the  sources  of  the  Incident  field,  produces  the  scattered 
field  (E^,H^)  exterior  to  S. 

In  the  absence  of  other  sources,  the  total  energy  radiated  by  the 
current  distribution  J  Is  always  positive  semldeflnlte: 

E(t)  -  J  W(t)  dt  >  0  ,  (A.2) 


where  E  Is  the  energy  and  W(t)  Is  the  total  power  radiated  by  J  at  time  t  and 
Is  computed  as  follows: 


W(t) 


(A.  3) 
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t^ere  the  lower  case  quantities  are  time-domain  counterparts  of  the  corre¬ 
sponding  upper  case  transform^ domain  quantities.  The  current  J  and  Incident 
field  e^  are  related  through  an  Integral  equation,  which  fact  will  be 
exploited  subsequently. 

We  next  suppose  the  Incident  electric  field  and  Induced  current  are 


given  by 


s^t 

s*t 

J(r,t)  »  3Q(r)e 

+  J*(r)e  "  -  2  Re' 

LV  “  ] 

S-t 

s*t 

r  s-ti 

e  (r,t)  -  ° 

+  Ej(r)e  -  2  Re 

(A. 4a) 


(A.  4b) 


where  s^  ®  >  0,  and  where  Eq  Is  an  arbitrary  complex  vector 

function  of  position,  and  Jq  Is  the  resulting  complex  current  response.  The 

excitation  Is  assumed  to  start  at  t  »  -®,  where  there  Is  no  Initial  energy 

<^0*=  -1 

In  the  system.  Since  e  =0  for  t  =  -«,  both  e  and  j  are  zero  and  there 
Is  no  transient  term.  That  Is,  the  'forced*  response  above  Is  the  total 
response. 

The  power  W(t)  radiated  Is  now 


and  hence  the  total  energy  radiated  Is 
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E(t)  ■Re 


f  fisLils.  +  °o  •  ^0 

■'s  1.  *0  “o  J 


If  we  write  the  first  Integral  In  polar  fom. 


S  "0 


then  E(t)  can  be  written  as 


/Vjo  if  Vi 

o  Srt  Ut. 


S  dS  e^"**  . 


S  “0 


(A.  5) 


J  -  dS  cos  (2u»-t  + 

s 


/®o  • 


>  0  . 


(A.  6) 


We  must  consider  the  case  ■  0  and  0  separately* 


Case  I,  Ug  ■  0 


In  this  case,  (A. 5)  can  be  written  as 


^  /  [Re  (Eg  .  Jg)  +  Re  (Eg  •  J*)]  dS 


'0  S 


2  ^ 
e  a«t 


r  /  (®o  •  «o  *  3;)J  is 


'0  s 


0  r 

—  I  Sil'o  •  "0} 


—  J  S.  E(,  .  R.  dS  , 
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from  which  we  conclude  that 


/ 


Re  E, 


Re  Jq  dS  >  0  . 


(A.7) 


But  since  Eq  Is  arbitrary,  we  could  replace  It  by  JEq  and  the  corresponding 
current  response  would  be  jJ^  by  the  linearity  of  the  system*  Hence,  from 
(A.7),  It  must  also  be  true  that 


/ 


Re  (JEq) 


RsCJJq)  dS  - 


Im  Jq  dS  >  0  . 


(A.  8) 


Equations  (A.7)  and  (A. 8)  together  Imply  that 

Re  /  Eq  •  j*  dS  -  /  (Re  Eq  •  Re  Jq  +  Im  Eq  •  Im  Jq)  dS  2  0  .  (A.9) 

S  S 

Case  II,  Uq  0 

From  (A. 6),  we  note  that  since  the  maximum  negative  value  for  the 
cosine  Is  -1,  we  must  have 


Re 


>  0  . 


(A. 10) 


Thus  for  both  cases  we  have 


s 


^0  >  0  . 


(A.  11) 


which  Is  a  necessary  condition  on  all  solutions  of  scattering  problems  involv¬ 
ing  passive  scatterers. 

Using  (A. 11)  we  next  derive  a  condition  similar  to  the  positive 
real  condition  on  driving-point  immlttances  In  circuit  theory.  The  condition 


applies  to  the  dyadic  kernel  of  the  operator  relating  J  to  E  through  the 
electric  field  Integral  equation 


n  x  E^(r,s)  ■  n  X  ^Z(r,r';s)  •  J(r';8)  ds'  . 


(A. 12) 


Since  a  and  J  are  transform  quantities,  (A. 11)  applies  for  a  point  a  *  Re  s 
In  the  right  half  plane: 


I  /j*(r)  •  Z(?,r')  •  J(?‘) 


dS  dS 


S  S 


h  J  f[j*(r)  •  Z(r,r')  •  3(r*)  +  [3*^)  •  Z(7,?’)  •  3(?’)]  dS 

S  s  •• 

%  /  •  Z(r,r’)  •  J(r’)  +  3*(r’)  •  Z+(?,?’)  •  J(^)]  dS’  d 

s  s 

h  J*J'[  J*(r)  *  Z(r,r’)  •  J<r')  +  J*(r)  •  Z*(r,r*)  »  J(r')j  dS'  dS 
S  S 

H  /  /[j*(r)  •  [ZU,^')  +  Z*U,1?')]  •  J(^’)]dS'  dS 

s  s 

//  J*(r)  •  Re  Z(r,r’)  •  J(r')  ds'  dS  >  0  ,  ( 


(A. 13) 


S  S 


where  the  dagger  denotes  the  transpose  conjugate  and  where  we  have  used  the 
reciprocity  condition  that 

zt(?,^’)  -  z*a,?’)  -  Z*(?’.?)  . 


Thus,  Re  Z  must  be  positive  semldeflnlte  for  a  >  0;  l.e. , 
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where  (A.  14)  Is  a  shorthand  notation  for  (A. 13),  except  for  operators  which 
are  merely  complex  constants. 

We  note  that  condition  (A.  14)  is  analogous  to  one  of  the  so-called 
positive  real  conditions  for  Immlttances  of  m-port  networks.  For  such  net¬ 
works,  the  positive  real  condition  is  shown  to  be  both  a  necessary  and 
sufficient  condition  for  the  realizability  of  a  passive  network.  The  suffi¬ 
ciency  of  the  condition  is  established  by  actually  deriving  algorithms  for 
synthesizing  such  networks.  We  are  not  yet  at  this  stage  in  electromagnetics, 
but  we  are  now  in  a  position  to  list  several  independent  conditions,  analogous 
to  those  of  an  m-port  immlttance  matrix,  which  the  impedance  operator  for  a 
passive  scatterer  must  necessarily  satisfy: 

(1)  Z(r,r';s)  is  analytic  in  Re(s)  >  0; 

(2)  Z(r,r';s)  la  real  for  real  positive  s; 

(3)  Re  Z(r,r';s)  >  0  for  Re  s  >  0. 

Though  the  sufficiency  condition  has  not  been  demonstrated,  we  nevertheless 
term  operators  with  kernels  satisfying  these  conditions  as  positive  real.  As 
in  network  theory,  it  is  possible  to  find  alternative  conditions  for  (l)-(3), 
which  are  simpler  to  test.  We  also  note  that  to  arrive  at  these  conditions, 
we  must  write  the  Integral  equation  in  the  form  of  (A. 12).  We  point  this  out 
because  it  is  comnxm  to  view  the  right-hand  side  of  (A. 12)  as  the  scattered 
field,  in  which  case  the  left-hand  side,  in  view  of  the  boundary  conditions, 
would  have  a  negative  sign.  In  this  form,  the  resulting  impedance  operator 
could  be  termed  'negative  real*. 
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Next,  we  note  thet  Z(r,r';s)  can  be  written  In  terms  of  its  eigen- 
spectrum  as 

Z(r,r';s)  -  ^  X^(8) J^(r;s) J^(r' ;s)  ,  (A.15) 


where  Is  an  eigenvalue  and  Is  the  corresponding  eigenvector  satisfying 
the  eigenvalue  equation 


/z(r,r';a)  •  J  (r*  ;s) 


dS  -  X^(a)J^(r;8) 


(A.16) 


and  where  Is  orthogonal  to  all  other  eigenvectors 


J  j_(r;s)  •  J  (r;s) 
«  m  n 


dS  »  6 


mn 


(A.  17) 


Note  that  Z  Is  an  operator-valued  analytic  function  of  s  and  we  assume  that 
X^  and  are  likewise  analytic.  The  absence  of  the  complex  conjugate  In  the 
second  eigenvector  In  (A. 17)  Is  at  first  disconcerting,  but  It  should  be 

ai 

remembered  that  the  operator  Z  Is  not  self-adjoint,  but  Is  complex  symmetric; 
l.e.,  Z+(r',r)  ■  Z*(r,r').  Briefly,  the  eigenvectors  should  be  blorthogo- 
nal  to  the  eigenvectors  of  the  adjoint  operator,  Z'*’(r',r),  which.  In  view  of 
the  complex  symmetry  of  Z,  are  just  J*.  The  resiiltlng  blorthogonallty  con¬ 
dition  Is  (A. 17). 

We  note  that  (A. 13)  must  be  true  for  all  possible  current  distribu¬ 
tions  J  and  hence  must  be  true  for  J(r;s)  »  J^(r;8).  Substituting  this  con¬ 
dition  and  (A.15)  Into  (A. 13)  results  In 
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Jj*(r)  •  Xx„(8)5  (r)J„(?’)  •  ds'  dS 

•;  o  *— *  n  n  n  m 

S  n 

-  Rej^X^Cs) 

-  Re^XJs) 

from  which  we  conclude  that  X^(s)  is  positive  real 

Re  X  (s)  i  0  ,  o  >  0  .  (A.  18) 

n 

It  can  be  shown  that  for  scalar  quantities,  (A. 18),  together  with  the  require¬ 
ment  X^(s)  Is  real  when  s  Is  real.  Is  equivalent  to  the  following  conditions: 

(a)  ^q(8)  has  no  poles  or  zeros  In  the  right  half  plane; 

(b)  poles  of  X^(s)  on  the  Imaginary  axis  must  be  simple  and  the  residues 
evaluated  at  these  poles  must  be  real  and  positive; 

(c)  Re  X^(j<i))  >0,  0<«<<». 


The  latter  condition  Is  particularly  important  because  it  requires  only  the 
examination  of  X^  on  the  ju  axis  rather  than  In  the  whole  of  the  right  half 
plane.  As  an  example  of  the  application  of  these  conditions,  we  note  that 
for  a  wire  loop  there  exists  only  a  4~dlrected  current,  which  Is  fotmd  by 
requiring  the  ^component  of  electric  field  to  vanish  on  the  loop.  Thus  the 
kernel  is  scalar-valued  and  has  been  found  by  T.  T.  Wu  to  be 
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Z(4,4')  -  I  jnoifCn  ® 


n— 


^  >]  ^  (A.  19) 

n-1 


where  the  rearrangement  to  the  trigonometric  form  is  made  to  more  closely 

compare  with  (A. 15).  Thus,  inust  be  positive  real.  But  ja^(s)  Is 

positive  real  If  and  only  If  is  positive  real  and  It  has  been 

1 


verified  by  direct  conputatlon  that 


jo^(s) 


has  no  poles  either  In  the 


right  half  plane  or  on  the  Jui  axis,  thus  satisfying  conditions  (a)  and  (b). 
Condition  (c)  implies  that  we  must  have 


Re 


Im 


Oj^(ju) 


i  0  , 


0  1  u>  S  “  , 


which  appears  to  be  true  [18]. 

Finally,  we  consider  the  ju  axis  poles  of  the  Inverse  operator, 
whose  kernel  Is  the  resolvent  kernel 


Z”^(r,r’ ,s) 


^(r,s)  jj^(r’  ,8) 


X,<s) 


(A.  20) 


We  know  that  poles  on  the  ju>  axis  must  be  simple  and  that  the  vector  compo¬ 
nents  of  the  current  which  are  cavity  mode  currents,  are  cophasal  and 
hence,  because  of  the  normalization  (A.17),  are  real.  Furthermore,  at  a 

pole,  one  of  the  eigenvalues,  say  X  ,  has  a  simple  zero  and  can  be  approxl- 

m 

mated  by  the  first  term  In  Its  Taylor  series  for  a  pole  at  Sq  by 
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X^<8)  -  Ct^(8  -  8q) 


or  in  polar  form  by 


X  (a)  -  a.e  , 
m  X 


where  (s  ~  8q)  ~  e  e^^.  Since  che  pole  term  will  dominate,  then  on  the 
semicircle  of  radius  e  in  the  right  plane,  we  have  the  condition 


Re 


/  fl*  •  Z  •  j 


dS  dS 


S  S 


If  j*(r,s)  •  J^(r,s)  J^(r',SQ)  •  JCr^.s^)  , 

C  C  M  ^ 


»  Re 


dS  dS 


S  S 


e 


/j  (r,s)  •  J(5,s)  dS 


2Re-~i 


Oj^e  e 


j* 


f J  (r,s)  •  J(r,8)  dS 


“ii' 


cos  (^  +  0)  i  0  , 


(A. 21) 


where  The  condition  (A. 21)  must  hold  for  -y  <  <  y  ,  l.e., 

for  8  in  the  right  half  plane,  which  is  possible  only  if  0  -  0.  This  implies 
that 

dX 


1  ds 


8“8, 


(A. 22) 


is  real  and  positive.  It  can  easily  be  checked  that  this  condition  is  con¬ 
tained  in  condition  (c)  above. 
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Note  that  for  anteana  problems .  we  can  take  to  be  a  unit  voltage 

source  (l.e. »  a  S  function)  at  some  pair  of  terminals  on  s  and  conclude  that 

the  input  admittance  is  positive  real.  The  above  proof  establishes  the  PRness 

of  the  inverse  eigenvalues,  but  does  not  address  the  question  of  the  PRness 

of  terminal  elgenadmlttances.  The  PRness  of  the  terminal  eigenadmittances 

may  be  explored  based  on  the  PRness  of  the  eigenadmittances  established 

n''*^ 

above  «ud  on  (2.13) ,  which  is  written  here  in  terms  of  the  natural  current 
modes  and  gap  geometry  as* 


-u 


d^<J  (r,.)  ;  J  (r»s)> 
n  n  ^ 


(A. 23) 


It  is  convenient  to  define  Yj^C®)  * 


F  (8) 


<J  (r,s)  ;  a  > 
a  g  S, 


4  <J  (r,s)  ;  J  (r,8)> 
n  n 


(A. 24) 


such  that 


i(8)  -  (yV  -  yM)  +  3(yV  +  yV)  . 


n  n  n  n 


n  n  n  n 


where  the  superscripts  denote  respective  real  and  imaginary  parts.  It  fol¬ 
lows  that  PRness  of  T^(s)  hinges  on  the  adherence  to 


yV  i  yV  .  Re{s}  i  0  . 


'n  n  'n  n 


(A.25) 


*The  foregoing  material  due  to  Wilton  used  his  notation  of  explicit  expression 
of  integration.  From  here  forward,  we  revert  to  the  syimetric  product  nota¬ 
tion  used  throui^ottt  the  body  of  this  work. 
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Because  of  the  complexity  of  (A. 24),  it  Is  difficult  to  draw  gen¬ 
eral  conclusions  regarding  the  satisfaction  of  (A. 25).  Clearly,  one  might 
test  (A. 25)  on  a  numerical  basis,  but  to  do  so  would  be  computationally 
costly  since  an  eigenvalue  problem  would  need  to  be  solved.  For  the  case 
that  the  elgenmodes  are  pure  real  for  s  ■  ju,  it  follows  directly  from 
(A. 24)  and  (A. 25)  that  the  terminal  elgenadmittances  Y^Cs)  are  PR.  This  is 
the  case  for  both  the  loop  and  sphere  geometries,  but  the  real  modes  for 
these  structures  devolve  from  symmetry  degeneracies.  On  the  other  hand,  the 
first  few  natural  modes  on  the  straight  wire  exhibit  small  real  parts  for  s 
near  the  ju  axis,  and  (A. 25)  is  likely  to  be  satisfied. 


With  the  respective  high-Q  and  low-Q  extremes  of  the  straight  wire 
and  the  sphere  likely  yielding  PR  terminal  elgenadmittances,  one  might  be 
tempted  to  draw  broad  conclusions.  However,  some  comaon  topological  feature 
such  as  convexity,  might  bear  on  the  results  for  these  two  special  cases, 
thereby  qualifying  any  general  conclusions  which  one  might  draw. 
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